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CHAPTER 1 
INTRODUCTION 


Smart composite structures offer the capability to combine the low density, superior 
mechanical and thermal properties of composite materials along with the inherent capabilities of 
smart materials to sense and adapt to their environments. Thus, the use of smart structures (also 
referred to as intelligent or adaptive structures) offers the potential to significantly improve the 
performance of aerospace structural components. However, before these materials can be 
implemented into actual structures, the coupled mechanical, electrical, and thermal behavior of smart 
materials must be fully characterized. This has led to extensive research since the 1980's to assess 
both the sensory and active responses of smart materials. 

1 . 1 Overview of Smart Structures 

Smart structures are distinguished from conventional structures by the presence of integrated 
actuator and sensor elements. In a typical smart structure application, the sensors are used to monitor 
the mechanical response of the structure through changes in the displacements, strains, or 
accelerations. Once an adverse or undesirable structural response is detected in the sensors, a 
controller generates the required input to the actuators. The actuators respond to this input and 
produce a corresponding change in the mechanical response of the structure to a more benign or 
acceptable state. The capability of smart structures to sense and adapt to their environment leads to 
a wide range of potential applications including: vibration suppression of aircraft structures; noise 
control of helicopter rotors; health monitoring of bridges; shape control of large space trusses; 
aeroelastic control of aircraft lifting components; and seismic control of buildings. Crawley (1993) 
and Loewy (1997) provide detailed overviews of the current state of smart structures research for 
aerospace applications. 

A variety of different materials can be utilized as either sensor or actuator elements in smart 
structure applications. Depending on the specific material used, the sensor and actuator elements are 
controlled through either electric, magnetic, thermal, or light energy. Some of the common actuator 
and sensor materials include: piezoelectric materials, shape memory alloys, fiber optics, 
electrostrictive materials, magnetostrictive materials, and electro-rheological fluids. Of the different 
materials available for use in smart structures, only piezoelectric materials have the unique capability 
to be used effectively as both actuator and sensor elements. Other advantages of piezoelectric 
materials which help account for their widespread popularity include: simple integration into the 
structure; a readily obtainable commercial supply of piezopolymers and piezoceramics; and 
familiarity in using these materials gained from previous applications in transducers. 

1.2 Piezoelectric Materials 


Historically, piezoelectric materials have been utilized mainly as active structural elements 
in transducers for application in strain gages, accelerometers, and sonar. Recently, the focus of 
research has shifted away from the transducer applications toward the development of smart structure 
applications which combine the active and sensory behavior of piezoelectric materials. 
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The basic characteristics of piezoelectric materials which allow for their use as sensors and 
actuators are the direct piezoelectric effect, converse piezoelectric effect, and the pyroelectric effect. 
In the direct piezoelectric effect, the application of a mechanical load on the piezoelectric material 
induces an electrical response. Through measurement of this electrical response, the mechanical state 
of deformation in the structure can be determined and monitored, leading to the sensory application. 
In contrast, the converse piezoelectric effect transforms an electrical input in the piezoelectric 
material into a corresponding mechanical strain. This leads to the active applications of piezoelectric 
materials, in which the state of deformation of the structure can be controlled or altered by applying 
the appropriate electrical input. The third characteristic behavior is the pyroelectric effect, in which 
the piezoelectric material responds to changes in temperature by producing an electrical response, 
which will influence both the direct and converse piezoelectric effects in changing temperature 
environments. 


1 .3 Objectives of This Research 

The renewed interest in the application of piezoelectric materials in smart structures has lead 
to extensive efforts to fully characterize the sensory and active behavior of these materials. Due to 
the complex coupled mechanical, electrical, and thermal behavior of piezoelectric materials, new 
analytical formulations and mechanics are required to accurately analyze the response of smart 
structures which contain piezoelectric elements. Thus, the purpose of this research is to develop 
comprehensive mechanics for accurately predicting the sensory and active behavior of composite 
laminates which contain piezoelectric elements. Special emphasis is placed on capturing the coupled 
response at the material level through the thermopiezoelectric constitutive equations and to introduce 
the displacements, electric potentials, and temperature as state variables in the analysis. This unified 
representation of the coupled response leads to the inherent capability to model both the sensory and 
active responses of piezoelectric composite laminates. 

The mechanics incorporate a layerwise laminate theory for more accurate analysis of 
displacements, strains, stresses, electric fields, and thermal fields, especially for thick laminates and 
laminates which contain strong inhomogeneties through-the-thickness. Thermal effects which arise 
from coefficient of thermal expansion mismatch, pyroelectric effects of the piezoelectric materials, 
and temperature dependent material properties are explicitly accounted for in the formulation. 
Corresponding finite element formulations are developed for piezoelectric beam, plate, and shell 
elements to provide a more generalized capability for the static and dynamic analysis of piezoelectric 
composite structures under arbitrary thermoelectromechanical loads and boundary conditions. The 
accuracy of the current formulation is verified by comparison with published experimental data and 
results from other analytical models. Additional numerical studies are also conducted to demonstrate 
additional capabilities of the formulation to represent the sensory and active behaviors. A future plan 
of experimental studies is also provided to characterize the high temperature dynamic response of 
piezoelectric composite materials. 


I 
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CHAPTER 2 

ANALYTICAL MODELING OF PIEZOELECTRIC MATERIALS 

2.1 Introduction 


The initial use of piezoelectric materials dates back to 1880, when the Curie brothers first 
discovered the direct piezoelectric effect. Until recently, the use of piezoelectric materials has been 
limited to a variety of transducer applications. The widespread use of piezoelectric materials as 
distributed actuators began only in the 1980's, when advances in design and manufacturing 
technologies made these applications feasible. The experimental work of Bailey and Hubbard (1985) 
is usually cited as the first application of piezoelectric materials as actuators for vibration control. 
Using a piezoelectric polymer film as the active element on a cantilevered beam, they were able to 
demonstrate active damping of the first vibrational mode. This new actuator application has led to 
renewed interest in the development of piezoelectric materials for advanced aerospace structures. 

As research into characterizing the active and sensory behavior of piezoelectric materials 
progressed, a variety of different analytical models were developed. These models can be classified 
into three broad categories as induced strain models, coupled electromechanical models, and coupled 
thermoelectromechanical models. The induced strain models use approximate theories to incorporate 
the piezoelectric effects and are generally limited to predicting only the active response of 
piezoelectric materials since the electric potential is neglected as a state variable in the formulation. 
The coupled electromechanical models provide a more consistent representation of both the sensory 
and active responses of piezoelectric materials by incorporating both the displacements and electric 
potential as state variables in the formulation. Typically, these models are implemented as finite 
element codes to provide a more general analysis tool and a wide variety of different beam, plate, 
shell, and solid elements have been developed. A natural extension of the coupled electromechanical 
models is to also incorporate thermal effects. These coupled thermoelectromechanical models 
include temperature as an additional state variable to account for thermal effects in addition to the 
piezoelectric effects. A more limited number of finite element codes have been developed with this 
capability. 


2.2 Induced Strain Models 


The induced strain models can be separated into two categories: (1 ) actuator models and (2) 
actuator and sensor models. The actuator models are concerned only with analyzing the active 
behavior of piezoelectric materials. They typically approximate the strain generated in the 
piezoelectric material by an applied electric voltage using statically equivalent forces and moments. 
The combined actuator and sensor models were developed to include predictions of the sensory 
response of piezoelectric materials. Although these models introduce the piezoelectric constitutive 
equations into their formulation, the electric potential is usually not included as a state variable, the 
conservation of electric flux is not considered in the equations of motion, and the sensory voltages 
are back calculated using the charge equation. 
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2.2.1 Actuator Models 


Crawley and de Luis (1987) developed an induced strain actuator model for beams. They 
formulated static and dynamic analytical models based on the governing equations for beams with 
attached and embedded piezoelectric actuators to model extension and bending. Experiments were 
performed on both isotropic and composite cantilevered beams with attached and embedded 
piezoelectric actuators to validate their models. The study found that segmented actuators are always 
more effective than continuous actuators since the output of each actuator can be individually 
controlled. They also showed that embedded actuators in composites degrade the ultimate tensile 
strength, but have no effect on the elastic modulus. Baz and Poh (1988) investigated methods to 
optimize the location of piezoelectric actuators on beams to minimize the vibration amplitudes. 
Numerical studies demonstrated the potential to control vibrations in large flexible structures using 
a small number of bonded piezoelectric actuators. Im and Atluri (1989) presented a more complete 
beam model which accounted for transverse and axial deformations in addition to extension and 
bending. Governing equations were formulated for a beam with bonded piezoelectric actuators for 
applications in dynamic motion control of large scale flexible space structures. 

Tzou and Gadre (1989) formulated an induced strain piezoelectric shell theory. A dynamic 
model was derived from Love's shell theory for application to multi-layered thin shells with active 
distributed actuators. A case study was validated with experimental results for the vibration 
suppression of a cantilevered beam with a piezopolymer actuator film. Crawley and Anderson 
(1990) developed a Bemoulli-Euler model to more accurately model actuation induced extension 
and bending in one-dimensional beams than the model of Crawley and de Luis (1987). The model 
neglected shear effects and was shown to be best suited for the analysis of thin beams and actuators. 
Clark et al. (1991) also extended the model of Crawley and de Luis (1987) to study the response of 
multiple piezoelectric actuators on beam excitations. Based on Euler beam theory, their model was 
validated with experimental results for the vibration response of a simply supported isotropic beam. 
The model was found to be best suited for performing initial studies to determine the optimal 
location of actuators for exciting specific vibration modes. 

Crawley and Lazarus (1991) developed induced strain actuation models for plates. Equations 
of strain actuation were derived for both isotropic and anisotropic plates. Exact solutions were found 
for simple geometries and boundary conditions, while approximate solutions were used to solve 
more complex problems. The models were verified with experimental results. Static analysis of a 
cantilevered composite plate with attached piezoceramic actuators were conducted to show the 
potential for shape control of structures. Dimitriadis et al. (1991) extended the one-dimensional 
induced strain beam model of Crawley and de Luis (1987) to two-dimensional plates with bonded 
piezoelectric actuators. Dynamic analyses were performed on simply supported plates to 
demonstrate the use of actuators to excite selective modes and the influence of actuator geometry 
on the modal response. Robbins and Reddy (1991) developed a piezoelectric layerwise laminate 
theory which was implemented into a beam element. Numerical comparisons were conducted using 
four different displacement theories (two equivalent single layer theories and two layerwise laminate 
theories) to demonstrate the increased accuracy in displacement and stress predictions obtained from 
using the layerwise theories. 

Wang and Rogers (1991) formulated an induced strain piezoelectric plate theory. The 
formulation was based on classical laminated plate theory and used Heaviside functions to represent 
the distributed piezoelectric actuators. The model was verified with the analytical work of 
Dimitriadis et al. (1991) and additional case studies were performed to demonstrate the capability 
of actuators to induce bending and extension in laminated plates. Clark et al. (1993) experimentally 
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validated the induced strain plate model developed by Dimitriadis et al. ( 1 99 1 ). Results from the free 
vibration analysis demonstrated capabilities to excite selective vibration modes of a simply supported 
plate using bonded actuators, as well as the significance of actuator location and excitation frequency 
on the structural response. Pai et al. (1993) accounted for geometric nonlinearities (large rotations 
and displacements) in piezoelectric composite plates. The developed partial differential equations 
characterize the dynamic response of composite plates through elastic and geometric couplings 
between the extension, bending, and twisting motions. 

Mitchell and Reddy ( 1 995b) formulated a power series solution for axisymmetric composite 
cylinders with either attached or embedded piezoelectric laminae. The solution was verified with 
finite element analysis. Numerical studies were performed to damp vibrations in truss-type 
structures using both an embedded cylindrical truss actuator element and an attached actuator patch. 
Lin et al. (1996) presented an induced actuation plate finite element based on first order shear 
deformation theory. Numerical studies were verified with analytical solutions and demonstrated 
capabilities to control the deflection of composite plates using piezoelectric actuators. Park and 
Chopra (1996) formulated one-dimensional models to predict the extension, bending, and torsion 
behavior of beams under piezoelectric actuation. Comparisons with experimental data showed good 
correlation only for applications in which actuators have low orientation angles (less than 45°) with 
respect to the beam neutral axis. 

Sonti and Jones ( 1996) developed differential equations of motion for a composite cylindrical 
shell with surface bonded piezoelectric elements. Approximate analytical solutions were obtained 
for the equivalent forces exerted by the actuator and shell. Chandrashekhara and Varadarajan (1997) 
implemented a finite element model for laminated composite beams with integrated piezoelectric 
actuators derived from a higher order shear deformation theory. Numerical studies investigated the 
effect of stacking sequence and boundary conditions on the actuator voltages and demonstrated the 
capability to achieve adaptive shape control of beam structures. Charette et al. (1997) formulated 
analytical models based on a variational approach to study plates with piezoelectric actuators. The 
model was specialized for a simply supported plate and used to determine the effects of piezoelectric 
actuator on the dynamic behavior of the plate. Results of the study were verified with experiments 
and showed only slight changes in the mode shapes of the plate. 

Chattopadhyay and Seeley (1997) implemented a third-order laminate theory into a finite 
element formulation to investigate piezoelectric actuators. Numerical results were verified with 
published experimental data. Additional comparisons were conducted to demonstrate the limitations 
of classical laminate theory in analyzing through-the-thickness stress and strain variations. Librescu 
et al. (1997) presented a model for composite beams with piezoelectric actuators that included 
structural tailoring and boundary-moment control. Numerical results demonstrated the potential to 
improve the dynamic responses of thin-walled cantilevered structures. Bhattacharya et al. (1998) 
developed a finite element formulation based on first order shear deformation theory for beams and 
plates. An eight noded isoparametric piezoelectric element was developed to perform free vibration 
analysis of laminated composite beams and plates. Numerical studies assessed the impact of stacking 
sequence, boundary conditions, and applied electric potentials on the free vibration response. 

Oguamanam et al. (1998) formulated a piezoelectric composite beam finite element using 
von Karman nonlinear strain-displacement relations. Numerical studies demonstrated the influence 
of stress stiffening effects on the natural frequency of slender beams. Tong et al. (1998) developed 
a two dimensional thin plate finite element to investigate shape control applications. Numerical 
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studies were performed to determine the optimum applied voltage, actuator layout, and actuator 
number for shape control of composite plates with distributed piezoelectric actuators. Chattopadhyay 
et al. (1999) implemented a third-order laminate theory into a finite element formulation to 
investigate the dynamic response of delaminated smart composite plates. Numerical studies were 
conducted on piezoelectric composite plates with single and multiple delaminations to demonstrate 
changes in the dynamic response. Hong and Chopra (1999) developed an induced strain plate finite 
element for composite plates based on classical laminate theory. The model was verified with 
experimental results and demonstrated the capability to achieve shape control using piezoelectric 
actuators. 


2.2.2 Actuator and Sensor Models 

Lee and Moon (1989) studied the control and sensing of bending and torsional deformations 
produced by an applied electric field using piezopolymer bimorphs. The experimental results 
validated the analytical model and demonstrated the use of piezopolymers as actuators. Lee (1990) 
developed a model that incorporated the piezoelectric constitutive relations. The model was based 
on classical laminated plate theory and was able to predict both the active and sensory behavior of 
piezoelectric materials. Lee et al. (1991) investigated the use of sensor actuator pairs for active 
damping control. Experimental studies were conducted on the active damping control of the first 
mode of a cantilevered plate using a sensor and actuator pair to validate the piezoelectric plate 
theory. 

Chandrashekhara and Agarwal (1993) developed a laminated piezoelectric plate element 
based on a first-order shear deformation theory applicable to both thin and moderately thick 
laminates. Numerical results were verified with previously published results for a cantilevered plate 
with attached piezoceramic actuators subjected to a static electric field. Hwang and Park (1993) 
developed a piezoelectric plate element based on classical laminate theory and Hamilton's principle. 
The piezoelectric constitutive relations were used to formulate a four noded, two-dimensional 
quadrilateral plate element for both sensory and active applications. Case studies were performed 
to investigate the static response of a piezoelectric bimorph beam and the vibration control of a 
cantilevered plate with attached piezoelectric sensors and actuators. Koconis et al. (1994a, 1994b) 
developed separate analytical models to investigate the sensory and active behavior of piezoelectric 
composite beams, plates, and shells. One model is used to predict the change in shape when a 
specified electric voltage is applied to the actuator, while the second model is used to determine the 
electric voltages necessary to achieve a desired shape. Both models were formulated using a two- 
dimensional, linear, shallow shell theory which includes transverse shear effects and validated with 
other numerical, analytical, and experimental results. 

Sung et al. (1996) derived sensor and actuator equations for a cylindrical piezoelectric 
composite shell. Based on classical laminate theory, these equations formed the basis of a sensor 
and actuator design methodology to control flexural and torsional vibrations in cylindrical shells. 
This methodology was used to design an experimental rig to demonstrate capabilities of the modal 
sensor and actuator to monitor and control the different vibration modes. Lam et al. (1997) 
developed a finite element model based on classical laminated plate theory for the active vibration 
control of composite plates with distributed piezoelectric actuators and sensors. Numerical studies 
were conducted on a cantilevered composite plate to demonstrate capabilities for static and dynamic 
analysis. Plettner and Abramovich ( 1 997) implemented a consistent methodology based on Kirchoff- 
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Love thin shell theory to model the static and dynamic response of anisotropic laminated 
piezoelectric shells. The formulation replaced the induced piezoelectric strain with an equivalent 
mechanical load. The model was verified with experimental and finite element results for a 
rectangular isotropic plate. 

Peng et al. ( 1998) implemented a beam finite element using a third order laminate theory for 
active vibration control of piezoelectric composite beams. Numerical studies were conducted to 
assess shape control applications and to investigate the effect of sensor and actuator locations on the 
response of the beam. Liu et al. ( 1 999) formulated a plate finite element based on classical laminated 
plate theory for modeling the static and dynamic response of composite plates containing 
piezoelectric actuators and sensors. Numerical studies were conducted to verify the model with 
results from previously developed models and to study the influence of stacking sequence and 
sensor/actuator position on the response of composite plates. 

2.3 Coupled Electromechanical Models 

Attempts to develop a more comprehensive representation of piezoelectric material behavior 
led to the development of more consistent models that captured the coupled response between the 
mechanical and electrical behavior. In these coupled models, the charge equation is incorporated into 
the equations of motion and the electric potential is introduced as an additional degree of freedom 
in the analysis. These models are generally also implemented as finite element programs to provide 
a more flexible and general purpose analytical tool. 

2.3.1 Analytical Models 

Mitchell and Reddy (1995a) formulated a refined hybrid theory for laminated piezoelectric 
composite plates. The displacement fields are modeled using third order shear deformation theory 
while electric potentials are represented using a layerwise laminate theory. An analytical solution 
was developed for simply supported boundary conditions and numerical results demonstrated the 
limitations of the induced strain methods in modeling thick laminates. Heyliger and Saravanos 
(1995) developed exact solutions to predict the vibration characteristics of simply supported 
laminated piezoelectric plates. Numerical studies were conducted to determine the influence of 
different laminations and aspect ratios on the natural frequencies and mode shapes. Batra and Liang 
(1997) developed three-dimensional elasticity solutions for the simply supported rectangular 
laminated plate with embedded piezoelectric layers. Numerical studies examined the steady state 
vibration of both thin and thick plates containing one actuator layer and one sensor layer. 

2.3.2 Finite Element Models 

The groundwork for much of the current research into finite element formulations was 
initiated by Allik and Hughes (1970) with the development of a three-dimensional tetrahedron 
element. Derived at a time when piezoelectric materials were used mainly as crystals in transducer 
applications, the finite element formulation incorporated the piezoelectric constitutive relations and 
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demonstrated the potential advantages of utilizing the finite element method. Naillon et al. (1983) 
formulated a finite element model to study the resonance phenomena of single piezoelectric 
structures typically used in ultrasonic transducers. Numerical studies were performed on the 
resonance characteristics of two-dimensional parallelepiped bars to assess potential applications in 
the design of ultrasonic probes. Lerch (1990) formulated two and three dimensional finite elements 
for performing vibrational analysis of piezoelectric sensors and actuators. The models were verified 
with experimental data. The natural frequencies and mode shapes of various piezoelectric structures 
were determined and used to optimize applications as transducers. 

Tzou and Tseng (1990) developed a thin piezoelectric solid element. Derived from 
Hamilton’s principle and the piezoelectric constitutive relations, the element is specifically 
formulated for thin plate and shell structures with distributed piezoelectric sensors and actuators. 
Numerical studies were performed on the vibration response of a cantilevered plate with both an 
active and sensory layer of polymeric piezoelectric material. Lammering (1991) developed a 
piezoelectric shell element. Based on the Reissner-Mindlin shell theory and incorporating the 
piezoelectric constitutive equations, a shell element was formulated for thin shell structures with 
attached piezoelectric layers. Case studies were conducted on a cantilevered beam with an attached 
piezoelectric polymer layer. Ha et al. (1992) formulated a three-dimensional brick element. Using 
a variational principle and the piezoelectric constitutive relations, they developed an eight-noded 
solid element. Results from static and dynamic case studies were verified with experimental results 
for composite plates with attached piezoceramic actuator and sensor patches. 

Heyliger et. al. (1994) implemented a layerwise laminate theory into a finite element 
formulation for plates. Two separate layerwise models were developed which incorporated the 
coupled equations of piezoelectricity to account for both the active and sensory behavior of 
laminated plates with piezoelectric layers. Numerical results for a simply supported composite plate 
with attached polymer piezoelectric layers were verified with exact solutions. Ray et al. (1994) 
developed a two-dimensional quadrilateral element using a higher order laminated plate theory. An 
eight-noded quadratic isoparametric quadrilateral element was formulated. Results from the static 
analysis of a simply supported cross-ply laminated plate bonded with a piezoelectric polymer were 
verified using previously reported exact solutions. Shieh (1994) developed a multiaxially active and 
sensory laminated piezoelectric beam element. Based on adjusted elementary beam assumptions to 
account for warping effects, the element can simultaneously model axial extension, biaxial bending, 
and torsional twisting of the beam. Numerical studies were performed on a space antenna frame to 
demonstrate capabilities for three-dimensional multiaxial vibration control. 

Saravanos and Heyliger ( 1 995) developed a layerwise finite element formulation for beams. 
Two separate theories were used to perform static and free vibration analysis of composite beams. 
Numerical results were verified with previously published analytical studies and demonstrated the 
increased accuracy of stress-strain predictions with the layerwise theory. Shen (1995) developed 
a finite element formulation for beams containing piezoelectric actuators and sensors. Based on 
Timoshenko beam theory, the methodology captured the coupling between the longitudinal and 
bending motions. The theory was validated with previously published analytical results for a 
piezoelectric polymer bimorph beam and experimental results for a cantilevered beam with attached 
piezoelectric sensors and actuators. Suleman and Venkayya (1995) formulated a plate element for 
analyzing composite plates with layered piezoelectric actuators and sensors. Based on classical 
laminate theory, a four-noded bilinear Mindlin plate element was developed. Previously published 
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experimental and analytical results for a cantilevered piezoelectric bimorph beam and a cantilevered 
composite plate with distributed piezoelectric actuators were used to validate the formulation. 

Donthireddy and Chandrashekhara (1996) also formulated a layerwise theory for beams. 
Results from the static response of a cantilevered composite beam with attached piezoelectric 
actuators was validated with previous analytical results. Additional parametric studies were 
conducted to study the influence of boundary conditions and ply orientation on the shape control of 
beams. Heyliger et al. (1996) implemented a layerwise laminate theory into a finite element 
formulation for shells. Results were verified with exact solutions for the static and free vibration 
response of a simply supported plate. Additional studies were performed on the active and sensory 
response of a cylindrical shell. Kim et al. (1996) developed a transition element to connect three- 
dimensional solid elements to flat shell elements. The reported approach used solid elements to 
provide detailed models of the piezoelectric material, while flat shell elements were used to provide 
a more flexible model for the substrate structure, and transition elements were used to connect the 
two regions. The model was verified experimentally for a cantilevered plate. 

Samanta et al. (1996) extended the eight-noded quadrilateral element developed by Ray et 
al.( 1994) for dynamic analysis. Based on a higher order shear deformable displacement theory, the 
model was developed for active vibration control of laminated plates with integrated piezoelectric 
layers. Numerical results were performed on a simply supported cross-ply plate with attached 
piezopolymer layers to demonstrate the potential to achieve significant reductions in vibration 
amplitude. Kim et al. (1997) provided additional details of the theoretical development of the 
transition element reported by Kim et al. (1996). Numerical studies were conducted to demonstrate 
convergence characteristics and to show the increased computational efficiency of this approach. 
Saravanos (1997) presented a shell element for curvilinear piezoelectric laminates which combined 
a first order shear deformation theory for the displacements along with a layerwise theory for the 
electric potential. The quadratic element was intended for static and dynamic analysis of thin to 
moderately thick shell structures. Numerical studies quantified the effects of curvature on the active 
and sensory response of piezoelectric shells. 

2.4 Coupled Thermoelectromechanical Models 

All of the previously described models neglect the implication of thermal effects on both the 
active and sensory response of piezoelectric structures. Although Mindlin (1974) derived the two- 
dimensional thermopiezoelectric equations for plates over twenty years ago, only limited research 
has been performed into this area. The development of models for thermopiezoelectric materials can 
be separated into two categories: (1) analytical models and (2) finite element models. The analytic 
models extend existing piezoelectric laminate theories to account for thermal effects and obtain 
solutions for specific problems. The finite element models provide a more general purpose tool to 
efficiently analyze complex problems. 

2.4. 1 Analytical Models 

Tauchert (1992) developed a thermopiezoelectric laminate plate theory. His theory extended 
classical laminate theory to account for thin laminated plates with thermopiezoelectric layers. 
Specific solutions were developed for both free and simply supported composite plates with attached 
piezoelectric polymer layers. Numerical results demonstrated the capability to reduce thermal 
deformations through application of active electric voltages. Tzou and Howard (1994) formulated 
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a thermopiezoelectric thin shell theory for applications to active structures. The generic shell theory 
was derived using Kirchoff-Love shell theory and Hamilton's principle. Using a simplification 
procedure based on the Lame parameters and radii of curvatures, specific solutions were obtained 
for a cylindrical piezoelectric ring, a piezoelectric ring, and a piezoelectric beam. Tang and Xu 
(1995) developed dynamic solutions for a simply supported anisotropic piezothermoelastic 
composite plate. The coupling between the elastic field and the electric and thermal fields were 
neglected to simplify the analysis. Numerical results demonstrated a significant reduction in 
deflection of a plate by the addition of a piezoelectric layer with a harmonic electric field. 

Tzou and Bao ( 1995) extended the thermopiezoelectric thin shell theory developed by Tzou 
and Howard (1994) for applications to anisotropic shell laminates with distributed sensors and 
actuators. The governing equations were simplified and applied to a thin piezothermoelastic 
laminated shell made of a piezoelectric polymer. Applications demonstrated the coupling between 
the elastic, electric, and thermal fields and the importance of all three fields on the overall behavior 
of the shell. Stam and Carmen (1996) presented axisymmetric thermoelectromechanical solutions 
for concentric piezoelectric cylinders. The analytical approach was used to model the quasistatic 
response of a linear piezoelectric motor. Results of the study demonstrated the capability to 
extrapolate the nonlinear dependence of the piezoelectric coefficients with electric fields to lower 
temperatures using the constitutive equations. Friswell et al. (1997) developed a linear model to 
investigate active damping of thermally induced vibrations. Numerical studies were conducted on 
a simply supported aluminum beam with a piezoelectric sensor/actuator pair to demonstrate the 
influence of pyroelectric effects on vibration control. 

2.4.2 Finite Element Models 

Rao and Sunar ( 1993) developed a finite element formulation with applications for integrated 
sensing and control of thermopiezoelectric materials. Numerical studies were performed on a 
piezoelectric bimorph beam and an isotropic beam with attached piezoelectric polymer layers. The 
results demonstrated the significance of thermal effects on the performance of distributed control 
systems. Jonnalagadda et al. (1994) implemented a nine-noded Lagrangian plate element using a 
first-order shear deformation theory. Numerical studies were performed on a simply supported 
composite plate with an attached piezoelectric active layer and demonstrated the importance of a 
higher order laminate theory to accurate predict shear deformations in thick laminates. Tzou and Ye 
(1994) extended the previously developed solid element of Tzou and Tseng (1990) to account for 
thermal effects. The resulting three-dimensional thin hexahedron element represented the 
displacements, electric potential, and temperature as state variables. Numerical studies were 
performed on a cantilevered isotropic beam with attached piezoelectric layers to demonstrate the 
influence of thermal effects on sensing and control. 

Chandrashekhara and Tenneti (1995) developed a nine-noded finite element for active 
thermal control of composite plates with piezoelectric actuators and sensors. The approach 
incorporated an induced strain approach to approximate the piezoelectric and thermal strains. 
Numerical studies demonstrated the capability to suppress thermally induced deformations through 
the application of electrical voltages in piezoelectric patches. Shen and Weng (1995) implemented 
a three-dimensional brick finite element to investigate composite plates with piezoelectric layers. 
Numerical studies demonstrated the significant coupling between the strain and electric fields and 
the capability to achieve thermal shape control of a simply supported piezoelectric composite plate. 
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Sunar and Rao (1997) formulated finite element equations for the design of thermopiezoelectric 
sensors and actuators. Numerical studies were conducted on cantilevered beams with piezoelectric 
actuators to show the significance of temperature effects on distributed control. 

2.5 Limitations of Existing Analytical Models 

The previously described analytical models contain three limitations that restrict the general 
application of these models. The first limitation is found in the induced strain models and arises from 
the use of approximate forces to represent the piezoelectric strains. This approximate representation 
fails to capture the coupled mechanical and electrical response and limits these models for use in 
predicting only the actuator behavior of piezoelectric materials. The induced strain limitations are 
overcome in the coupled electromechanical and thermoelectromechanical models through the use 
of a more consistent representation of the coupling which occurs between the electrical and 
mechanical responses. However, both the coupled electromechanical and thermoelectromechanical 
models are restricted by a second limitation arising from the predominant use of single layer laminate 
theories to predict the laminate response. In general, the single layer theories only provide good 
predictions for thin laminates and are inaccurate for both thick laminates and laminates which 
contain strong inhomogeneities through-the-thickness. The third limitation is found in both the 
induced strain and coupled electromechanical models, and arises due to neglecting thermal effects, 
even though these effects can be significant in changing temperature environments. Even in the 
coupled thermoelectromechanical models which account for thermal effects, many of the analytical 
formulations fail to account for all the mechanisms which produce thermal effects in piezoelectric 
materials: induction of thermal strains due to thermal expansion mismatch, pyroelectric effects, and 
temperature dependent material properties. Currently, no analytical model addresses all three of these 
limitations. 

Thus, in this research, each of these three limitations are addressed to provide an enhanced 
analysis capability for modeling arbitrary structures which contain piezoelectric materials. A 
consistent representation of the coupled mechanical, electrical, and thermal response of piezoelectric 
materials is implemented to overcome the induced strain limitations and provides an inherent 
capability to model both the sensory and active behaviors of piezoelectric materials. A layerwise 
laminate theory is implemented to address the single layer laminate theory limitations and leads to 
more accurate analysis of thick laminates and laminates with contain strong inhomogeneities 
through-the-thickness. Thermal effects arising from thermal expansion mismatch, pyroelectric 
effects, and temperature dependent material properties are explicitly accounted for in the 
formulation. The development of this enhanced formulation has been reported for a linear beam 
element (Lee and Saravanos, 1996), a bi-linear plate element (Lee and Saravanos, 1997, 1998) and 
a quadratic curvilinear shell element (Lee and Saravanos, 2000). 
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CHAPTER 3 

THEORETICAL DEVELOPMENT 


This section outlines the foundations and steps for developing the analytical formulation for 
piezoelectric composite materials. This development is based on the theory of linear piezoelectricity 
which have been previously reported in detail by Cady (1964), Tiersten (1969), Nye (1972), and 
Parton and Kudryavstev (1988). 

3.1 Current Formulation Assumptions 

The current formulation is based on the linear theory of piezoelectricity in which the 
equations of linear elasticity are coupled to the charge equations of electrostatics through the 
piezoelectric constants. Due to the assumption of a linear theory, both mechanical and electric body 
forces and couples are neglected in the derivation. The current model also assumes infinitesimal 
deformations and strains, and, thus, no distinction is made between the initial and final positions. 
The current formulation also neglects the two way thermopiezoelectric coupling and focuses only 
on the one way heat conduction analysis. This assumption essentially provides that temperature 
changes due to changes in strains and electric fields are small compared to the magnitude of the 
thermal load. Thus, temperature changes only produce mechanical and electrical forcing on the 
piezoelectric material. 


3.2 Governing Material Equations 

The mechanical response of the piezoelectric material is represented by the equation of 
motion, 

°ij.j + fi = P ", ( 1 ) 

while the electrical response is described by the electrostatic equation for the conservation of electric 
flux, 

D u = 0 ( 2 ) 


where Oy are the stresses, fj are the body forces per unit volume, p is the density, u, are the 
displacements, and D ; are electrical fluxes. In addition, i, j ranges from 1 to 3, superscript dots 
represent time derivatives, and subscript commas represent differentiation. Through use of the 
divergence theorem and neglecting body forces. Equations (1) and (2) can be expressed in an 
equivalent variational form as 


/ (p 


5 u. 


a.. 

v 


S V 


D 5 £.) dV = f t i Sir dr + J q 5 <p dT 


( 3 ) 


where Sy represent the strains, E ; are the electric fields, t, are the surface tractions applied on the 


N AS A/TM — 200 1 -2 1 0892 


13 



surface T t , q is the electrical charge applied on the surface T p of the piezoelectric material, cp is the 
electric potential, and V represents the whole volume including both composite and piezoelectric 
materials. 


3.1.1. Strains and Electric Field in Cartesian Coordinates 

Two additional relationships for the strain and electric field are required in the equivalent 
variational form (Equation 3). In a Cartesian coordinate system, the small deformation strain- 
displacement relations are 

S = — (u. . + li .) (A\ 

ij 2 J ' 1 K 

while the electric field vector is related to the electric potential by 

E i = ' %. (5) 

3.1.2. Strains and Electric Field in Curvilinear Coordinates 

The curvilinear coordinate system employed in the formulation is shown in Figure 1 . Each 
ply of the laminate remains parallel to a reference curvilinear surface A 0 and an orthogonal 
curvilinear coordinate system is defined, such that the axes £ and rj lie on the curvilinear 
reference surface A 0 , while the axis £ remains straight and perpendicular to the layers of the 
laminate. A global Cartesian coordinate system O xyz is used to define A c , hence, a point r = (x,y,z) 
on the curvilinear laminate is, 

r(^,Ti,0 = r 0 (^ri) + ^ (6) 

where, r 0 = (x 0 ,y o ,z o ) are the Cartesian coordinates of the reference surface A 0 and £ indicates the 
unit vector perpendicular to the reference surface. 

The strain and electric field in this curvilinear coordinate system take the following form 
(Fung, 1965) 

s „ - -7= [<W> - «*< r ' i * r ‘'>] (7) 

2 V S >i g jj 


( 8 ) 
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Figure 1 : Curvilinear piezoelectric laminate and coordinate systems. 


where g u are the metric tensors and Ty are the Christofel symbols. The following form of the 
Euclidian metric tensors are employed in the current formulation (Soedel, 1993), 

Sn = (1+ 4“) 2 *n (9) 

R \ 

*22 = O *22 (10) 

I\^ 


*33 * 


(ID 


where R, are the local radii of curvature and g n ° are the components of the metric tensor on the 
surface A„ ( C=0) defined as 
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Substituting Equations (9)-( 1 1 ) into Equation (7 ) produces the following relationships (Soedel, 1 993) 


'ii 


O O 

1 , 8 u. n 8 ii , 

(u , + -V-’ + VV’) 

o., 9 x " 9 ° 

Sll(l + -r-) 822 

n. 


(14) 


1 


o o 

<?22,£ §22 

( U + V’ + VV ) 

r o -I /? 

S \ v. . 


gl1 


(15) 


5 33 = H ’.4 


(16) 


c 1 , #22.c \ ^ 1 

■ 7 — 7 V > * _ 

S ,-; ( |.A) S„ SnO*-£-) S22 


, ^ll.n x 

— (V,- -«) 

s 






(17) 


5 13 = M .c + 




o 

— (w - - 


(18) 


S„ = v, + 


1 


*22< 1+ T-) 


O 

(vv - — — v) 

C X ' n ^ 


(19) 


where u, v, w are the displacements of the curvilinear system. Similarly, introducing Equations (9)- 
(11) into Equation (8) produces the following relationships 




*4 


(20) 
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(21) 


E, = 


822(1 + i~) 


— ‘P.n 


E 3 = -9 C 


( 22 ) 


3.3 Thermopiezoelectric Constitutive Equations 

The constitutive equations for a linear thermopiezoelectric material employing standard 
contracted notation (Nye, 1964) are 

s. = 4 ' T <o “n * £ „ ♦ “Tin e 

(23) 

- 4L(T) <y„ * C r (0 E t * p° m r m e 

(24) 

or in semi-inverted form, 


= Cg(T) S, -*LiT)E m - 6 

(25) 

D,„ = eJjT) S„ * 4 r <r) £, + p s /(T) 9 

(26) 

with the following relationships in the various piezoelectric and thermal properties arising from the 
inversion, 

el<T) = cgcn <<r) E m ( 27 ) 

X e /(T) = C^(T) a E /(T) 

(28) 

g(T) = *£<r) - eJ(D 

(29) 

/’"(D = P,/(T) ~ d E a (T) X E /(T) 

(30) 


where C a(J and s ft(! are the elastic stiffness and compliance tensors; d ([ni and e (mi are the different forms 
of the piezoelectric tensor due to inversion; e in); is the electric permittivity tensor; a u and X a are the 
different forms of the coefficient of thermal expansion due to inversion; p m is the pyroelectric 
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constant; 0 = AT = T - T 0 is the temperature difference from the current temperature T and the 
thermally stress free reference temperature T 0 ; and; superscripts E, a, S, and T, represent constant 
voltage, constant stress, constant strain, and constant temperature conditions, respectively; a, p = 1, 
...,6; and k, nr = 1,2, 3. 


3.3.1 Piezoelectric Beam 

The specific form of the thermopiezoelectric constitutive equations for a piezoelectric beam 
are simplified since only axial variations of the displacements are assumed. 


p. ] 


O 
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0 d 3] 
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E 3 . 
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(31) 


^3 *31 a i + e 33 ^3 + Pi ® 


(32) 


or in semi-inverted form 
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0 c 
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0 € 


31 


e i5 ® 
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(33) 


^3 e 3l ^1 + £ 33 + ^3 ® 


(34) 


3.3.2 Piezoelectric Plate and Shell 

For the piezoelectric plate and shell, the thermopiezoelectric constitutive equations for a 
monoclinic orthotropic piezoelectric material (Mason, 1950) are used. 
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(35) 
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CHAPTER 4 

LAMINATE MECHANICS 


In order to provide more general and accurate analysis of piezoelectric laminates, different 
theories are implemented for the beam, plate, and shell laminates. A layerwise laminate theory is 
used as the basis for the different formulations since this theory represents a generalization of other 
single and higher order laminate theories (Reddy, 1 993). The layerwise theory allows separate fields 
to be assumed for each state variable which provides for more accurate analysis of thick laminates 
and laminates containing strong inhomogeneities. In the current formulation various simplifying 
assumptions and approximations are used for the beam, plate, and shell laminates. 

4.1 Layerwise Laminate Theory 

The general layerwise laminate theory is formulated by introducing piecewise continuous 
approximations through-the-thickness for each state variable. 


N 

u(x,y,z, t) = J^u j (x,y,t) \\i J u (z) (39) 

Li 

AT 

v(x,y,z,t) = X v j (x,y,t) ^ J y (z) (40) 

7-i 

N 

w(x,y,z, t) X w J (x,y,t) \| f£,(z) (41) 

j - 1 
N 

(p (x,y,z,t) = X^U'^O V f (z) (42) 

7 = 1 
N 

Q(x,y,z,t) = X V 0 (s) (43) 

7 = 1 


where N is the number of interpolation functions y and u J , v j , w j , (p 1 , and O' are the generalized 
laminate state variables. In general, any order of interpolation functions can be used and different 
interpolation functions can be selected for each state variable depending on the level of detail desired 
in the analysis. 


4.1.1. Beam Laminate Theory 

In deriving the beam laminate theory, only axial and through-the-thickness variations of the 
state variables are assumed. In addition, the displacement w is assumed to be constant through-the- 
thickness. This leads to the following simplified approximations from the general layerwise theory. 
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( 44 ) 


N 

U (X,z,t) = Y l > J i- x ^) V ( " ) 

j - 1 


w(x,z,t) = w°(x,t ) (45) 

N 

9 (x,z,t) = Y V(x>0 Viz) (46) 

7 = 1 
N 

d(z,r) = Y ^(O V(2) (47) 

7 = 1 

¥ , U) = [1-^1 (48) 


where N is the number of interpolation functions i|/ j , £=z/h, and h represents the thickness of each 
discrete layer. In this formulation, the same linear Lagrangian interpolation functions are used for 
each state variable. 


4.1.2 Plate Laminate Theory 

The plate laminate theory is formulated with the displacement w assumed to be constant 
through-the-thickness, which leads to the following simplified approximations from the general 
layerwise theory. 


N 


u(x,y,z,t) = Y u Viz) 

7=1 

(49) 

N 

v(x,y,z, t) = Y v J (x,y,t) Vi z ) 

7 = 1 

(50) 

w(x,y,z,t) = w°{x,yj) 

(51) 

N 

(p(jc,y,-,r) = Y <p / (*»;y>0 Viz) 

(52) 


7 = 1 
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( 53 ) 


N 

Q(ZJ) = E P'O) V(2) 

j i 

V(z) = [l-5 5] < 54 > 

where N is the number of interpolation functions V )/ 1 and c=z/h. For the current plate laminate theory, 
the same linear Lagrangian interpolation functions are used for each state variable. 

4. 1 .3 Shell Laminate Theory 

The shell laminate theory combines linear displacement fields through-the-thickness of the 
laminate for the displacements u and v (along the £, and rj axes respectively) with layerwise electric 
potential and temperature fields through the laminate, consisting of N discrete continuous segments. 
Consequently, the present shell theory is specialized for modeling the deformation of thin and 
moderately thick piezoelectric shells, while maintaining the capability to capture the through-the- 
thickness electric and thermal inhomogeneities. The various field variables are approximated by the 
following form. 


«<5,n&f)= « 0 ( 5 ,v) + 5^(5, v) 

(55) 

v&il&')= v’&V) + 5P„(5.V) 

(56) 

w(4,r|,^,t)= w°(£,iy) 

(57) 

<p(5,tl,^) = XV&v) ^(Q 

j- 1 

(58) 

e(£,o,L,n- E 

j - i 

(59) 

V(0 = [1-5 51 

(60) 


where N is the number of interpolation functions \p J , ^=Cfh, and (L and P n are the flexural rotation 
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angles. The same linear Lagrangian interpolation functions are used for the electric potential 
and temperature state variables. 

4.2 Generalized Laminate Matrices 


The next step in the development of the formulation is the introduction of the layerwise 
approximations into the variational form of the equation of motion and the electrostatic equation 
(Equation 3). By performing an integration through-the-thickness, the dependence of the through- 
the-thickness coordinate can be separated from the in-plane coordinates into generalized laminate 
matrices for the density, stiffness, piezoelectric, dielectric permittivity, thermal expansion, and pyro- 
electric properties. 


4.2. 1 Beam Laminate Matrices 

The equivalent variational form is derived by substituting Equations (4), (5), (20), (21), and 
(31)-(34) into Equation (3), and integrating through-the-thickness. 


^2 f Pn H m t>u k dx + f p 33 w 5 w dx + 

A-’l m-1 ' 


IE/( D* 

k 1 m - 1 J 


km du m dbu k r\ km ms: k \ j 

+ £> ss m 5m ) dx + 

dx dx 55 


E C n k , d vv ^ ; d 8 w , , 

/ B 55 ( Sm + u k ) dx + 

k - 1 J dx 


dx 


r . 5h> 55w , 

I a *j; — dx * 


dx 


ee r £,‘r iv" 

k - 1 m-1 J 


55m k du m s k \ j 
+ 5(p ) dx - 

dx dx 


E E f Gg r 5(p*' dx £ f (A* ^ - P k 5<p* ) dx 

k - 1 m I J k 1 J OX 

x x 

= ft. duJr + fg 5 (pdr 

r r p 


(61) 
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in which the dependence of the z-coordinate has been separated into the generalized density matrix 
[p], the stiffness matrices [A], [B], [D], the piezoelectric matrix [E], the dielectric permittivity matrix 
[G], the thermal expansion matrix [AJ, and the pyroelectric matrix [P 3 ]. 

The laminate density matrices are 


Pn = £ f b P V'"(-) dz 

/_1 J Z 


L. 

P 31 = £ f b p dz 

/ l J Z 


the laminate stiffness matrices are 


A S5 f b C 55 dz 

l 1 J ~ 


B 


55 


E / b c 55 * 

l - 1 J: dz 


D n = £ f b c uV*(z) v|/ m (z) rfz 

oi" -£ *£& dz 

i=\ J z dz dz 

The piezoelectric laminate matrix is 

r, km r , v 5v)/ w (z) * 

E 31 = b b e 3l V (2) ~T~ * 

/-I J c <?Z 

The dielectric permittivity laminate matrix is 

< = £ f : h % ^ ^ dz 


(62) 


(63) 


(64) 


(65) 


( 66 ) 


(67) 


( 68 ) 


(69) 
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The laminate thermal force vector is 


L 

A , = E ( b\ & V(~) vj/* ( 3 ) dz 
/ 1 J 


The laminate thermal electric displacement vector is 

= E f ^ yu) * 

/ = 1 J OZ 

x 

where L is the number of plies and b is the width of the beam. 
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4.2.2. Plate Laminate Matrices 


The equivalent variational form is derived by substituting Equations (4), (5), (24)-(25), and 
(36)-(40) into Equation (3), and integrating through-the-thickness 


N N 


E E f (Pn + P 22 v m Sv k )dA + f p 33 *8wdA 

k = 1 m--\ \ \ 


N N 


E E f + ♦ V.;vl) + 

k = 1 m \ \ 
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A km / m o k m s k m k \ a . . m s: , . ^ . 

A 16 (w v 8 m + M v OM v + li x 0 V v + l’ v 0 M V ) + A 22 V y 0 V v + 

x km , trt k m o A' w r* k m ^ \ 

A 26 (m v 5v v + V v 8w v + v r 8v v + V, 5v v ) + 


a£(u”&u* + m.;6v; + v a m 8v>; + v“5«;) }<m 


jV 


E f {BZ(v m 8W' V + H’ V 8v w ) + 5«(n m 5w v + + 
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E E f ids 

k = 1 m = l 


km v m bv k + D,**(h m 5v * + v m 5« *) + D™u m &u k + 


. kfti / ms;,. * 


«*™ m s,, A 


£ 3 T«.rv' + £ 3 T(v.;v + ^(<v + v;v> + 

£ 3 Y'(p m 5«j' + EgqTfo* + £3^(9 m 8w.y + <P m 5vJ) 

- GS"q>;&p* - } dA 
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Eft- A78«; - A™ 8 v * + P, m S< + P 2 " 8 q»; + ^3 8 <p m } dA 

m = 1 \ 

A 

= J f. 8w ( .dT + | q 5<prfT 


(72) 


in which the dependence of the z-coordinate has been separated into the generalized density matrix 
[p], the stiffness matrices [A], [B], [C], [D], the piezoelectric matrix [E], the dielectric permittivity 
matrix [G], the thermal expansion matrix [A], and the pyroelectric matrix [P]. 
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The laminate density matrix is 


P f = E f P ¥"'(c) dz 

ii J 


(73) 


for i j = 1 1 and 22, while 


L „ 

p 3 3 = E ( p dz 

1 1 j 


(74) 


the laminate stiffness matrices are 


4 - 


L , 

E / 

/- i J 


c,. y*(z) ¥'"(") 


(75) 


for ij = 1 1, 12, 16, 22, 26, 44, and 66; 


B,r - 1 / C„ * 

/=i J J dz 


(76) 


C 

y 


L. 

E f 

l l J 


C,j * 


(77) 


D 


r = e / 

/=i J 


c d\|/*(z) dv/^z) 
^ dz dz 


(78) 


for ij = 44, 45, and 55. The laminate piezoelectric matrices are 

L 


e? = E 

1-1 J 1 dz 


(79) 


= E f ru) r<z) d 

1 1 J 1 dz 


(80) 
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for ij = 31, 32, and 36. The laminate dielectric permittivity matrix is 


Gjj"' = £ f e y*(z) Viz) 

i - i J 


* 


( 81 ) 


for ij = 1 1 and 22, while 


r km _ r ( Z ) cV" ( C ) , 

^33 2^ I e 33 3 3 

/ - 1 ' OZ OZ 


(82) 


The laminate thermal force vector is 


A? = £ f \ Vj °j /(-) * 

/ i • 


(83) 


for i = 1 and 2. The laminate thermal electric displacement vector is 


P- = £ f P, V/ 0 y Y*(z) * 

M J 


(84) 


for i =1 and 2, while 




(85) 


where L is the number of composite plies and piezoelectric layers. 

4.2.3 Shell Laminate Matrices 

The equivalent variational form is derived by substituting Equations (11), (12), (24)-(25), and 
(40)-(45) into Equation (3), and integrating through-the-thickness 
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J { ( p ' it " + p s p,) da ° + (p- 4 v° + p B (3 n ) Sv ° + p‘ 4 vr "8vv "+ 


(P B H° + P°P 4 ) 5p, + (p B v° + p y P n ) 5P n }^ n + 




N 


E / 

./‘I J . 


S" + A 12 S 2 " + + 5 12 A- 2 + - A\Gt)&S° 


(A 12 5," + A 22 5 2 ° + B yi k { + 5 22 * 2 + £ 3 V - A 2 6 / )55 2 ° + 

J 

(^4 + 4— >*4 + <A*4 ^ B 66 k 6 )8S 6 ° 7 - (5 66 5 6 " + D 66 k 6 m 6 } 


8 11 


(5, ,5" + B i2 S 2 ° + £>„*, + Z> 12 * 2 + 4^’ - Ai0 / )8*i }<Mn 
(fi 12 S," + £ 22 S/ + D n k { + D 22 * 2 + 4^ “ A 2 0 , ')5^ 2 + 


* 


A' f N C ° rn A/ 

E E f {(4— - 4 — ^ + />/*— )5<p* + 

j~ [ kl j > e,, (£,,) g n 


.4 


O 

8 1! 


ten) 


(4— - G. 


5/ 


jk ‘P.n n jk O' „ * 

+ Pi )8(p„ 

n ' 1 


£22 (£ 11 )” £22 

(44 + 44 + 4*1 + 4^2 - 4/ + 4* 0t ) 5 «p / }^rfn 

= j t j 5w ( .</r + J y 5(pjr 


( 86 ) 


in which the dependence of the z-coordinate has been separated into the generalized density matrix 
[p], the stiffness matrices [A], [B], [C], [D], the piezoelectric matrix [E], the dielectric permittivity 
matrix [G], the thermal expansion matrix [A], and the pyroelectric matrix [P]. 

The laminate density matrix is 

P^ = £n £22 E/P^ (87) 


P B = £h£22 E ftP<% 


( 88 ) 


P° = £u£22 E jV P^ 

c 


(89) 
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The laminate stiffness matrices are 


=*,"*£ E f c ^ 

l 1 J 


(90) 


for ij = 11, 12, 22, 44, 55,66 


B,, - SuSn E f 

l 1 J . 


( 91 ) 


0„ = Sn«22E /C 2 ^ 

M •(. 


(92) 


for ij = 11, 12, and 22. The laminate piezoelectric matrices are 


e; = E f % **«;>«*; 

/ --i •( 


(93) 


for ij = 14, 15, and 24. 


= <?n'g 2 2 E f e ij 

i-\ i 


VCO 

X, 




(94) 


= Sn822 E f 
/ = 1 , 


<V(Q 




(95) 


for ij = 31, 32, and 36. The laminate dielectric permittivity matrix is 

c,‘" = SuSn t f h V*<0 V"(0 <<C ( 96 ) 


for ij = 1 1 and 22, while 




V(0 arco d r 

x, % 


( 97 ) 
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The laminate thermal force vector is 


A? = 8 u 822 E f \ d(, 

M J 

L 

Af - Snfe E [ V*(0 dl 

i - 1 J , 


for i = 1 and 2. The laminate thermal electric displacement vector is 

P/"' = SnS 2 2 E f P k V*(0 

/ - l 



L 

o o C 

S \\&22 2 s / 
M 


/>3 


a; 


for i =1 and 2, while where L is the number of composite plies and piezoelectric layers. 


(98) 


(99) 


( 100 ) 


( 101 ) 
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CHAPTER 5 

FINITE ELEMENT FORMULATION 


Finite element based solutions for structural problems are obtained by incorporating addition- 
al local in-plane approximations as detailed in Hughes (1987) and Cook et al. (1989). For 
piezoelectric materials, the following approximations are made for the generalized laminate state 
variables introduced in Equations (39)-(42), 


M 


'(x,y, t) = N‘ u (x,y) 


i 1 


M 


v j (x,y,t) = Y v ji (t) N‘Xx,y) 


(102) 


(103) 


M 


M 


VV 


yj 


'(x,y,t) = Y w Ji (t) N' v (x,y) 


i i 


(104) 


M 


<p '(x,y,t) = Y^(t) K(x,y) 


i 1 


(105) 


M 


&(x,y, t) = Y N ^(x,y) 


i-i 


(106) 


where M is the number of in-plane shape functions N. Typically, any order of shape functions can 
be used and different shape functions can be selected for each state variable, depending on the level 
of detail desired in the analysis. 

By implementing the in-plane approximations for the state variables into the variational form 
of Equation (3) and collecting the coefficients, the following discrete matrix form can be obtained. 


[M ] 0 

L UU J 

{«} 

> 1 

[*j (Kj 

( <«)L 

W)} 

0 0 

. {<p}j 

T 

wy [kj. 

1 { 9 }] 



[KJieil 

i^eliei) 


(107) 


where the submatrices K uu , and indicate the elastic, piezoelectric and permittivity matrices; 
K^e and are the thermal expansion and pyroelectric matrices of the structure; M uu is the mass 
matrix; F are the applied mechanical loads; and Q are the applied voltages. 

The coupled finite element formulation can also be expressed in a compact form with the 
electric potential partitioned into active and sensory components such that 
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( 108 ) 


[M ] 

L UU J 

0 

({*}) 


tO 

f{«) 1 

0 

0 


K ] 

tO 

W}J 


inoHOl^l-KJie} 


W(p 


rFA 




tp(H 


where superscripts F and A indicate the partitioned submatrices in accordance with the sensory (free) 
and active electric potential components, respectively. Thus, the left-hand side includes the unknown 
electromechanical responses of the structure, u and <|> F (i.e. the resultant displacements and voltage 
at the sensors, respectively). The right-hand includes the known excitations of the structure by the 
mechanical loads, applied voltages on the actuators, applied temperature loads, and electric charges. 

The matrix equations can also be condensated into the following independent equations for 
the sensory electric potentials, 

i<p f )=[0 l <[O{^vtv{ 0 >-[C F ]iw-{e''M)) do?) 


and the structural displacements, 

WJ{U] +([AJ -[<]lO '[ClH") = 

{ F(t ) } - ( [<f] [O 1 O + [O ) { ^ ) + (110) 

[< F : ] [< F i l {Q F (n}-( [<f ] [ O 1 [k^ + [kj ) { e j 

5.1 Beam Element 

The piezoelectric beam element formulation is obtained by incorporating the following in- 
plane approximations to the state variables, 

M 

u j (x,t) = £ u ji {t) N‘(x) (111) 

1 = 1 

M 

w j (x, t) = £ w ji (t) N'(x) (11 2) 

1=1 

M 

•p'uo = IV’W (113) 

i=i 
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N'(.x) = [i-q g 


(114) 


where M is the number of in-plane linear shape functions N, ^=x/L, and L is the length of the 
element. A selectively reduced integration scheme is also implemented for the second stiffness term 
(containing D„) of Equation (67) in order to eliminate locking. 

Incorporating Equations (111 )-( 1 1 4) into Equation (61) results in the following specific 
forms of the submatrices defined in Equation (107) for the beam element, where the submatrices are 
defined using the beam laminate matrices of Equations (62)-(71). The mass matrix is 


= l l A 'l r 


[p„i 

0 



[N] dx 


(115) 


The stiffness matrices are 


[*J - f x 


[A y r [£„] [Af,J + [N] t [D 55 ][N] [^] r [5 55 ] [N, x ] 


[N] t [B 55 ][N, x ] 




dx (11 6) 


IV = 

f [N] T [E 3l ][N, x ] dx 

(117) 

IV 

= f [A] r [A,] dx 

J X 

(118) 

IV = 

f [N] T [G 33 ][N] dx 

(119) 

IV 

= f [N] t [P 3 ] dx 

(120) 
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5.2 Plate Element 


The piezoelectric plate element implements the following in-plane approximations to the 
state variables. 


M 


u j (x,y,t) = w 71 ( t) N'(x,y ) 


1 1 


( 121 ) 


M 




( x,y,t ) = Y v Ji (t) N'(x,y) 


M 


( 122 ) 


M 

w J (x,y, t) = £w 7 '(Y) N'(x,y) 

i - 1 


(123) 


M 

<pi(x,y,t) = N‘(x,y) 

1-1 


(124) 


/V'u,y) = I [(i- 4 )(i-n) (i + s)(i ~ri) (i+^)(i+n) (i-^)(i+n)] (125) 

4 


where M is the number of in-plane bi-linear shape functions N , £=x/a, and r|=y/b. 

Implementing Equations ( 1 2 1 )-( 1 25 ) into Equation (72) leads to the following finite element 
submatrices for the plate element, where the submatrices are defined using the plate laminate 
matrices in Equations (73)-(85). The mass submatrices are 


wj = f A w 


tp„i 

0 

0 


0 

[ P221 
0 


0 

0 


[N] dA 


[P 33 ] 


(126) 


The stiffness submatrices are 


[AT ] = 

L uu 1 


[*i,J 

[^ 12 ] 

[^, 3 ] 


[^1 2 ] 

1^ 22 ] 

[K^] 


[^ 3 ] 

[^ 2 3] 

[^33] 


(127) 


where 


A 


[N, x ] t [A u ][N, x ] 

[N] t [D 55 ][N] 


+ [MJ^[A 16 ][7V, v ] + + 

- [N. y ] T [A"\[NJ } dA 


(128) 
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[K n ] 


f { [N, x l T [A l 2 ][N. y ] + [Ag r [A 16 ][M,] * [M v ] r [A 26 ][iV, v ] + 

A 

[N] T [D 45 ][N] + [iV, v ] r [A 66 ][M v ] } dA 


[AT 13 ] = / { [N] T [B 45 ][N, y ] + [M r [fi 55 ][MJ } 

A 

[tf 22 ] = f { [M V ] T [A 22 ][M V ] + [JV,/[\][^] + [M v ] r [A 26 ][yv, v ] + 

A 

[NV[D m )[N] * [M v ] r [A 66 ][A^, v ] } dA 


[K a ] = f { [A'] r [5 44 ][M v ] + [Aq r [fi 45 ][yV,J } dA 

A 


[AT 33 ] = f { [M/tCJIAg + [M v ] r [C 45 ][M v ] 

A 

[MJ r [C 55 ][MJ } dA 


lKf[C 45 ]tN, x ] + 



[/V, v f[£ 31 ][iV] 

[Ay r [£ 32 ][tf] 


[A^, v ] r [£ 36 ][ATl 

[MJ r [£ 36 ][An 



[MJ r [A 2 ] 


( 129 ) 


( 130 ) 


( 131 ) 


( 132 ) 


( 133 ) 


( 134 ) 


( 135 ) 
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(136) 


[* w ] = / ( [M,] r [G n J[M v ] + [M v ] r [G 22 ][A5 v ] + [yV] r [G„][A^] } dA 


[*„,] = f { [M,] r [/ > 1 ] + [N,f[P 2 ] + [Af]V 3 ] } dA 


5.3 Shell Element 


The piezoelectric shell element implements the following in-plane approximations to the 
state variables, 


M 

Uj°(^r],t)= J^Uj°\t)N'(& i) 


/=] 


(138) 


M 

P„(5.V)= EP«(G(V'(5,n) 


;1 


(139) 


M 

<p w (5.v)= E<p' m (G (v'(5ji) 
/ = 1 


(140) 


M 

e m (5,v)= E0 w '(G(v'(^t d 

/-i 


(141) 




(1-5)0 -n) d-5 2 )(i-n) - G-5 )(i-ti 2 ) 

(i + 5)d-n) - (i-5 2 )d-n) - 0 + 5Xi-n 2 ) 

(1+5X1+TD - (i+5)(i-n 2 ) - (i-5 2 )0+n) 

(i-5)(i + t )) - (i-5 2 )(i + n) - (i-5)(i V) 

(i-5 2 )d-n) 

(l+5)(l-ri 2 ) 

(1-5 2 )(1+ti) 

(1-5)(1-T1 2 ) 


(142) 


where M is the number of in-plane shape functions N for an eight-noded Serendipity element. 
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By introducing Equations (138)-( 142) into Equation (86), the following finite element sub- 
matrices for the shell element are calculated in accordance with the shell laminate matrices defined 
in Equations (87)-( 101 ). The mass submatrix is 


[MJ - I [AH 3 

J A 


[p^l [Pg] 

[P B ] [P D 1 


[N] dldx\ 


(143) 


The stiffness submatrices are 


where 


[tf ] = 

*- LIU J 


[*.l] 

[* 12 ] 

[^13] 

[*J 

[* 1 5 ] 

[ K n ] 

[* 22 ] 

[tf 23 ] 

[*24] 

[* 24 ] 

[^13] 

[tf 23 ] 

[^33] 

[J*] 

[* 35 ] 

[*J 

[^4] 

[^34] 

[J*] 

[* 4 5 ] 

[^13] 

[^25] 

[^35] 

[*« ] 

[*55] 


(144) 


[ K u ] = f { - J —[ N , t Y [ A n m ,] 

{ (8i i) 2 

° \2 


TrA hm 1 . 2g22 ’£ [^] r [A 12 ][M £ ] + 


(8n 822 ) 2 


[An r [A„][An 




/ 0 \2 ° 
( 8 n) 822 


[A^f[A 55 ][(V] 


(822) 


— ([N,J-^[N]) t [A 66 ]([N,J-^[N]) } d^dr\ 
22 ) 2 8n 8n 


(145) 


N AS A/TM — 200 1-210892 


39 


1 


+ 


o 

§ 11 ’ 


[tfp] = f { — ■ 1 ^-[N„] T [A U ][N] + -L- 

{ (S,y-g22 8 M 822 


[/V,,] r [A I 2 ][M n ] 


g||,n ^ [(Vf[A n ][A(] + 


O 

8 22 > 


- n 


0 o v ^> 


(§ II §22 )■ 

1 


« / 0 \2 
8 11 (§22) 


[W] r [A 22 ][W, n ] + 


O 0 

. aN ] ([ N ]-^[N]) } dUr] 

O O ■ O S 0 

&11 <?22 &11 <?22 


[^„] = f { — [ttgfl A U ][N] + — ^— [M e f[A p ][^] 
a §,1 R ] ' 8 n * 2 


<?22’£ 


o o D 

§11 §22 R i 


■ [yvf[A n ][/v] 


o 

8 22 ’c 


o 6> r> 
8 11 <?22 ^2 


■[Wfl^HW] 


5 11 8 22 


•([M n ] -§n, n [AH) r [A 55 ][^] } d^d n 


9 0 

[K l4 l = f { — - — [M £ ] r [fi n ][iV, £ ] + [M £ ] r [g, 2 ][A^] + 

a (§11 > 2 (§ ii ) 2 §22 ' 

- ( * - 4> [Aq r [g„][Ar] - ^-[^f[A 55 ][^] + 

(§11 §22) 1 

o o 

— ([N, n ]-^[Nlf[B 66 ]([N ]-^[N]) } d^dx] 

O 1 O ” o 

822 §11 §11 


( 146 ) 


(147) 


(148) 
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[K„] = / ( 


o 

8 ii, 


/ o \2 ° 

(§n) §22 


'^-[N,,] T [B n ][N] + — — [Ay r [* 12 ][M n ] + 


O O 

<§ 11,11 <§ 22.6 

/ o o . 2 
(§11 § 22 > 


[^] r [ 5 n ][^] 


O 

8 22., 


0 z 0 x 2 

§11 (§ 22 ) 


0 0 

§11 §22 


[ An r [ fl 22 ][ M n ] + 


— [N]) } d^d n 

0 0 1 0 n 0 

§11 §22 § 1 1 §22 


/ 0 \2 

[*r„i - f ( d*™L 

— J , 0 o>,2 

A (§11 § 22 ^ 


[N] T [A..][N] + -^^[N] T [A n ][N,J 


§n (§22> 2 


12 J L' ' ’ll J 


1 


< 0\2 
(§ 22 ) 


[MJ r [A„][AU + -2-[A^] 7 '[A 44 ][Ar] + 


(*,r 


0 a 

— ([M £ ]- — [A(]) r [A 66 ]([MJ- — [AH) } 


(§ 1 .) 


o 
§ 22 


§22 


° a ° 

[IK*] = f { -^-[N] T [A n ][N] + ‘ ^L -[An r [A I2 ][An - 

‘ J A §, 1 § 2 ° 2 *, §11 §22 *2 

— — [M n ] r [A 12 ][/V] + [M n ] r [A 22 ][A^] - 

§22 ^1 §22 ^2 

— — [A^A*] [tf] } d<^df\ 

§22 *2 


(149) 


( 150 ) 


( 151 ) 
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( 152 ) 


l* 24 ] = / { 


O 

§11 


’n 


o o 


[A(] r [fi u ][A r , £ ] + gll,n g224 [Nf[B n ][N] + 


.4 (5ll)'5 2 2 


(5,, *22 ) 


1 


0 0 'I 

5ll 522 


[M B f[fl l2 ][M £ ] 


0 

§22. c 


o , 0 v 2 

5] 1 (522 1 


— ([(V, £ ]-^[(V]) T [5 66 ](W ]-2i2H[Ar|) ) <%d n 

0 0 ^ 0 wu ’I o 

5 1 1 522 522 5 1 1 


[M„n* 22 ][An + 


5ii’ 


[^25] = f i --- l n> [iV] r [g,,][AT] 


2 5 , M 


j ✓ 0 

4 (§11 § 22 ) 


0 z 0^2 
§11 ( §22 ) 




— 7TtM n f[5 22 ][iV, n ] - ^[(Vf[A 44 ][An + 
(§22 ) 


R, 


0 0 

([M e ]- — [(V]) r [S 66 ]([)V, c ]-^[)V]) } dtdr] 

y 0 0 v 0 " 0 

( 5l 1 522 1 522 522 


(153) 


[^34] = f {-^fUnlM + r- 4 -[(V] r [A 12 ][iV] + 

A R ' R 2 

- -4r[M n ] r [A 44 ][M n ] + 
r: (5 22 ) 

— — [M»] r [A 5 ,][MJ } <%<h 1 

,0,? *= 

(5 U )- 


(154) 
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1 


o 

IK U ] - f { [JV] r [B„][Ay - - [W] r [B| 3 ][W] 

J . <*,”> R ' *° R. 


&II 8 2 


— — [Nf[5 p ][iV, t ] + -JB±—[M] t [B 22 ][N) + 
8 11 ^2 #1 1 S22 R 2 

-^[M/[A 55 ][An } d\d n 

Sn 


+ 


0 

[*„] = f { glKn ~ [N] t [B u ][N] * -l-[N] T [B l 2 ][N, n ] + 
‘ { «,**»*, * 22*1 

— ^ — m 7 [fi P ][iV] + — — [^] r [fi„][M n ] + 

* 22*2 ‘ * 22*2 

-^[M n ] r U 44 ][An } d^dr] 

822 


[K u ] = f {[yvf[A, 5 ][iV] H- —L ~[N,J T [D U ][N,J + 

A (Sll)" 

2g2 ?i -[jV f[D 12 ][AG + ( — [^V] r [£> 22 ][AG - 

( 8 n ) 2 822 (Su S22) 2 

— ([^J-5n n [An) r [£> 66 ]([M n ]-^n. n [An) ) d\dx\ 
(S22) 2 


(155) 


( 156 ) 


( 157 ) 
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J 


o 

§ 11 


[AT 4S J = f {— 2 H*—[N„] T [D n ][W] 
J /„ °\2 „o 


J / l '\S l/ ' O o 

A (§[l) §22 §11 §22 


[H,J T [D l 2 ][N, n ] + 


O 

l,il S 22 .C rA/iT'rn u \n , ^22,c , *, .* ■ # » , ■ n , ■ + 


, o 0.2 

(s 11 <?22^ 


[Ny[D { 2 )[N] 


0 / 0 \2 
8 11 (<? 22 ) 


[N] r [D^][IV] 


o o 

£l] £22 


g ll 'lrAnx7-rn 1 / r aa i - ^ [N]) } d$dt\ 


IN])' [D 66 ]([N,J 


§11 


o 

§22 


/ ° \2 

( 8 1 1 . 1 ] ) 


[K„] = [ im T [A 44 ][N] + — L i J - [N] T [D n ][N] 

J / 0 °\2 


/ 0 0 \2 
(£ll £22^ 


2g,l - n [MJ r [£> 12 ][(V] 


0 \2 


(§ 22 > 


§U (§ 22> 2 " 

— aN , n ]-^[ N ]) T [ D 66 ]([N ] -^5 

() > () uu 1 O 

£ll £22 £22 


[ AU r [ D„][AU 


[ W ]) } d ^ di ] 


[K ^=L 


r [^ 3 |] T § 22 , t t ^ 32 ^ 

[ M £ ] r — ^[ AT ] + [^] r ’ 


o 

§11 


o o 
§11 §22 


T [ £ isl 

m - m T — — w 

O rj 
§11 


.8nJE 3i ] 


- vn , ^ " 3 l J T [£^1 T [ Eja ] 

[N] T = [N] + [N,J T -^-[N] - [N] T ^^[N] 

o S -, O ‘I 


§1. 


§22 


§22 *2 


[AT( 


[£■ 31 ] _ . r»r \T ^ ^24 ] r An , r 


R, 


•)[ n]+[n y 

Ix-f 


° \2 


(§ 22 ) 


-[N] + [N, t V-^-[N] 

- / °\2 

(§n) 


r [£ ls ] r [£ 31 ] 

[N] t —^~ [TV] + [W,J — [TV] + [TV] 7 32 [jV] 


o 

§11 


o 

§11 


§11 §22 


O O O T I 


§22 


§11 §22 


o 

§22 


dcdl] 


( 158 ) 


( 159 ) 
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[K„] = / 1 


[*,(]=/ (W 


rA.] g 2 2„[A 2 ] 

[Wf -[AU + [N] T i [N] 


O 

S n 


o o 

8 1| §22 


[A] r gU '" [Al] [N] + [MJ r — [M 


o o 

Sll S22 


822 


[AT(- 


[A, ] [A,] 
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CHAPTER 6 

CASE STUDIES AND DISCUSSION OF RESULTS 


Results from representative problems for a variety of different beam, plate, and shell 
structures are presented in this section. The accuracy of the current analytical formulation is verified 
with comparisons from published experimental data and other analytical models. Additional 
numerical studies are also conducted to demonstrate the capabilities of the current formulation to 
model the sensory and active behavior of piezoelectric composite materials. 

6.1 Comparisons with Room Temperature Experimental Results 

Results from previously published room temperature experimental studies are used to verify 
the developed analytical formulation. The problems examined consist of a cantilevered piezoelectric 
bimorph beam examined by Lee and Moon (1989) and a clamped composite plate with attached 
discrete piezoelectric patches studied by Crawley and Lazarus (1991). 

6.1.1 Piezoelectric Bimorph Beam 

The cantilevered piezoelectric bimorph beam problem examined by Lee and Moon (1989) 
is shown in Figure 2. The bimorph is created by bonding two layers of a piezoelectric polymer 
together such that the poling directions are opposite in each layer. In this configuration, the 
application of an external voltage will result in a pure bending deformation. The beam has a length 
(L) of 8 cm, a width of 1 cm, and a height of 220 pm. The material properties of the polyvinylidene 
fluoride (PVDF) piezoelectric polymer used in the experiment are listed in Table 1 . A finite element 
analysis of the bimorph beam problem is performed using the developed beam, plate, and shell 
elements using sixteen elements along the length of the beam with two discrete layers through-the- 
thickness. All three elements produced similar results, so only one set of predictions are presented. 
Figure 3 illustrates the predicted free end (at x/L = 1) displacement obtained by applying active 
voltages of between 0-500 Volts, along with the experimental results of Lee and Moon ( 1 989 ). Good 
agreement is observed between the predicted results of the current analytical formulation and the 
experimental data. 


6. 1 .2 Composite Plate with Piezoelectric Patches 

The cantilevered composite plate problem examined by Crawley and Lazarus (1991) is 
shown in Figure 4. The plate is constructed using AS4/3501 Graphite/Epoxy with a length (L) of 

29.2 cm, a width (b) of 15.2 cm, and a thickness of 0.83 mm. Two laminate configurations, [0/±45], 
and [30/0 2 ] s , are examined with fifteen G-l 195 piezoceramic patches with a thickness of 0.25 mm 
bonded on both the top and bottom surfaces of the plate. Figure 4 also depicts the finite element 
mesh used in the analysis: a 16 by 9 mesh with 8 discrete layers (one for each composite and 
piezoelectric layer through-the-thickness) is used for both the plate and shell elements. Only one set 
of results is presented since similar results are obtained using each element. The material properties 
of the composite and piezoelectric material are listed in Table 2. The purpose of the experimental 


N AS A/TM— 200 1 -2 10892 


47 



studies was to investigate the resulting deflections of the plate when active voltages are applied to 
the piezoceramic patches. The experimental displacements along the centerline (y/b=0.5 ) of the plate 
are shown in Figure 5 for the [0/±45] s plate with an applied voltage of 394 V/mm and in Figure 6 
for the [30/0 2 ] s plate with an applied voltage of 472 V/mm. Good agreement is observed between the 
predictions of the current analytical formulation and the experimental results. 



Figure 2: Cantilevered piezoelectric PVDF bimorph beam (Lee and Moon, 1989). 
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Table 1 : Material properties of PVDF piezoelectric polymer (Lee and Moon, 1989 ). 


Property 

PVDF 

Elastic Modulus, E (10 9 Pa): 

2.0 

Poisson's Ratio, u: 

0.29 

Shear Modulus, G (10 9 Pa): 

0.8 

Piezoelectric Charge Constants (10 12 m/V): 

23.0 

d 32 

3.0 


5 |— - - 

~ Current Formulation 


£ 4 i- * Lee & Moon (1989) 

a / 



0 100 200 300 400 500 


Electric Voltage (volts) 


Figure 3: Comparison of tip deflections of a PVDF bimorph beam under applied 
active electric voltages of 0-500 Volts at room temperature. 
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y 



Piezoceramic Graphite/Epoxy 0.83 mm 



~29.2 cm 


Figure 4: Cantilevered AS4/3501 composite plate with 30 attached discrete 
G-l 195 piezoelectric patches (Crawley and Lazarus, 1991). 
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Table 2: Material properties of piezoceramic (G-l 195) and 
Graphite/Epoxy (AS4/3501 ) (Crawley and Lazarus, 1991). 


Property 

Piezoceramic 

Graphite/Epoxy 

Elastic Moduli ( 10 9 Pa): 

E„ 

63.0 

143.0 

E„ 

63.0 

9.7 

E„ 

63.0 

9.7 

Poisson's Ratio: 

1) 12 

0.3 

0.3 

1^23 

0.3 

0.3 

d 3 i 

0.3 

0.3 

Shear Moduli ( 10 9 Pa): 

G p 

24.2 

6.0 

G 73 

24.2 

2.5 

G 31 

24.2 

6.0 

Piezoelectric Charge Constant (10 12 

m/V) : 


^31 

254. 

— 

^32 

254. 

— 

Electric Permittivity (10 9 f/m): 

£1 i 

15.3 

— 


15.3 

— 

£ 33 

15.0 

.... 
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Centerline Displacement (w/b) 


0.060 

Current Formulation 

0.048 [- 

• Crawley & Lazarus (1991) 


0.036 b 

0.024 

0.012 

0.000 - 

0.00 0.20 0.40 0.60 0.80 



f 



J 


1.00 


Axial Length (x/L) 


Figure 5: Comparisons of centerline displacement of an AS4/3501 [0/±45J s plate with 
attached G-l 195 piezoelectric patches under 394 V/mm at room temperature. 
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Figure 6: Comparisons of centerline displacement of an AS4/3501 [30/0 2 ] s plate with 
attached G-l 195 piezoelectric patches under 472 V/ram at room temperature. 
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6.2 Comparisons with Analytical Results for Temperature Effects 


Results from previously published high temperature analytical studies are used to verify the 
developed analytical formulation. The problems examined consist of a cantilevered piezoelectric 
bimorph beam examined by Sunar and Rao (1997) and a simply supported composite plate with 
attached discrete piezoelectric patches investigated by Ha et al. ( 1 992). These problems present slight 
variations on the room temperature experimental studies of Lee and Moon (1989) and Crawley and 
Lazarus (1991) to analytically quantify the influence of temperature effects. 

6.2.1 Piezoelectric Bimorph Beam 

The cantilevered piezoelectric bimorph beam problem examined by Sunar and Rao (1997) 
is shown in Figure 7. The bimorph has a length (L) of 100 mm, a width of 5 mm, and a height of 1 
mm. The material properties of the polyvinylidene fluoride (PVDF) piezoelectric polymer at 27 °C 
and at 77 °C are listed in Table 3. An applied electric voltage of 250 Volts was applied to the 
bimorph beam at 27 °C to induce a bending deformation. The impact of thermal effects on the 
deflection of the bimorph is subsequently examined by applying a uniform 50°C temperature 
increase. The current finite element analysis is performed using 20 elements along the length of the 
beam with two discrete layers through-the-thickness. All three elements (i.e. beam, plate, and shell) 
produced similar results, so only one set of predictions are presented. Figure 8 shows the predicted 
deflection of the bimorph along the length of the beam under an applied active voltages of 250 Volts 
at 27 °C and 77 °C, along with the corresponding analytical results of Sunar and Rao (1997). Good 
agreement is observed between the two sets of results, which demonstrate the impact of temperature 
effects on the deflection of the bimorph beam. 

6.2.2 Composite Plate with Piezoelectric Patches 

The simply supported composite plate problem examined by Ha et al. (1992) is shown in 
Figure 9. The [0/±45] s plate has a length (L) of 37.2 cm, a width (b) of 22.8 cm, and a thickness of 
0.75 mm. Fifteen equally sized G- 1 1 95 piezoceramic patches with a length of 6 cm, a width of 6 cm, 
and a thickness of 0. 1 3 mm are bonded on both surfaces of the plate. Figure 9 also depicts the finite 
element mesh used in the analysis: a 16 by 10 mesh with 8 discrete layers (one for each composite 
and piezoelectric layer through-the-thickness) used for both the plate and shell elements. Only one 
set of results is presented since similar results are obtained using each element. The material 
properties of the composite and piezoelectric material are listed in Table 4. Analytical studies were 
performed to investigate the capability to achieve thermal distortion management by applying active 
voltages to the piezoelectric patches to eliminate the bending deflection induced by applying a 
100°C linear thermal gradient through-the-thickness of the plate (50°C on the top surface and -50°C 
on the bottom surface ). The predicted displacements along the centerline (y/b=0.5) of the plate along 
with corresponding applied active voltages to produce this deflection are shown in Figure 10, along 
with the corresponding results from Ha et al. (1992). Good agreement is observed between the two 
predictions. The results demonstrate the capability to achieve shape control of the thermally induced 
bending deformation by the application of increasing active voltages to the piezoelectric actuators. 
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Figure 7: Cantilevered piezoelectric PVDF bimorph beam (Sunar and Rao, 1997). 
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Table 3: Material properties of PVDF piezoelectric polymer 
(Sunarand Rao, 1997). 


Property 

27 °C 

77°C 

Elastic Compliance, c n (10 9 Pa): 

3.8 

1.8 

Thermal Expansion, a (10 6 m/m°C): 

150. 

215. 

Piezoelectric Constant, e^ (C/m 2 ): 

0.046 

0.049 

Electric Permittivity, e ( 10 9 f/m): 

0.1026 

0.1553 

Pyroelectric Constant, P (10 2 C/m 2o C): 

0.04 

0.055 

Reference Temperature, T„ (°C): 

27.0 

27.0 


0.25 
| 0.20 
I 0.15 

£ o.io 

Q 

.& 0.05 

H 


Current Formulation 
Sunar and Rao ( 1997) # 


77 °C 






27 °C 


0.00 

0 2 4 6 8 10 


Axial Distance (cm) 


Figure 8: Comparison of tip deflections of a PVDF bimorph beam under an applied 
active electric voltages of 250 V at 27 °C and 77 °C. 
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Piezoceramic Graphite/Epoxy 0.75 mm 



Figure 9: Simply supported AS4/3501 composite plate with 30 attached discrete 
G-l 195 piezoelectric patches (Ha et al., 1992). 
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Table 4: Material properties of piezoceramic (G-l 195) and 
Graphite/Epoxy (AS4/3501) (Ha et al„ 1992). 


Property 

Piezoceramic 

Graphite/Epoxy 

Elastic Moduli (10 9 Pa): 

Eh 

63.0 

150.0 

E„ 

63.0 

9.0 

E„ 

63.0 

9.0 

Poisson’s Ratio: 

«12 

0.3 

0.3 


0.3 

0.3 

*6l 

0.3 

0.3 

Shear Moduli (10 9 Pa): 

G P 

24.2 

7.1 

G 2 3 

24.2 

2.5 

g 3 . 

24.2 

7.1 

Piezoelectric Charge Constant (10 12 

^31 

m/V) : 

254. 



254. 

— 

Electric Permittivity ( 10 9 f/m): 

Ell 

15.3 

— 

&>-> 

15.3 

— 

e 33 

15.0 

— 

Thermal Expansion (10' 6 m/m°C): 

Oil 

0.9 

1.1 

a 22 

0.9 

25.2 

«33 

0.9 

25.2 

Reference Temperature, T 0 (°C): 

20.0 

20.0 
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Figure 10: Comparisons of centerline displacement of an AS4/3501 [0/±45] s plate with 
attached G-l 195 piezoelectric patches under a 100°C thermal gradient through-the- 
thickness (50°C on top, -50°C on bottom) with 0, 31, and 61 V applied active voltages. 


NASA/TM — 200 1-210892 


59 



6.3 Comparisons with a Commercial Finite Element Analysis Program 


Evaluations of the developed analytical formulation for both the sensory and active responses 
are verified with corresponding results from a commercial finite element analysis program, 
ABAQUS (1996). In order to represent the location of piezoelectric layers, the standard laminate 
notation is expanded, such that the piezoelectric layer is indicated by the letter p. The problem 
examined consists of a Carbon/Epoxy composite plate with an attached layer of piezoceramic 
material. The [0 8 /p] plate has a length (a) of 50.8 mm, a width (b) of 25.4mm, and a height (h) of 
5.08 mm, as shown in Figure 1 1. The plate is clamped on one side and is subjected to a uniform 
thermal load of 50°C. The material properties used are listed in Table 5. A 20 x 10 mesh with two 
discrete layers (one for the Carbon/Epoxy and one for the piezoelectric layer) is used for both the 
plate and shell elements from the current formulation, while a 10 x 5 x 2 mesh is used for the 20 
noded continuum elements in ABAQUS. Only one set of results is presented for the current 
formulation, since similar results were obtained for both the plate and shell elements. Pyroelectric 
effects are intentionally neglected (i.e. set to zero) in the current formulation, since ABAQUS 
neglects such effects. 


6.3.1 Sensory Mode 

The piezoelectric layer can be utilized in either a sensory or active mode. In the sensory 
mode, the piezoelectric layer is free to develop a corresponding electric potential in response to the 
thermally induced deformation of the plate through the direct piezoelectric effect. The resulting 
displacements along the centerline of the plate (along y/b = 0.5) are shown in Figure 12, while the 
corresponding electric potentials which develop in the piezoelectric layer are shown in Figure 13. 
There is good overall agreement between the current formulation and the ABAQUS results. The 
differences in the sensory electric potentials arise due to the consideration of the through-the- 
thickness Poisson's effect in ABAQUS, which are neglected in the current formulation. In typical 
health monitoring applications, these sensory voltages would be used to provide the essential 
feedback necessary to infer the current state of deformation in the structure. 

6.3.2 Active Mode 

The piezoelectric layer can also be operated as an actuator by applying a voltage differential 
to the piezoelectric layer and utilizing the converse piezoelectric effect to produce a corresponding 
structural deformation. Figure 14 shows the centerline deflection (along y/b = 0.5) for three different 
applied electric potentials (0, 250, and 500 Volts). Good agreement is observed between the current 
formulation and the ABAQUS predictions. The results also demonstrate the potential of piezoelectric 
actuators to completely alter the initial thermally induced deformation of the plate by the application 
of increasing electric potentials in order to achieve thermal shape control. 
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Table 5: Material properties of piezoceramic and Carbon/Epoxy. 


Property 

Piezoceramic 

Carbon/Epoxy 
(Daniel and Ishai, 1994) 

Elastic Moduli (10 9 Pa): 

E„ 

69.0 

142.0 

E^ 

69.0 

10.3 

e m 

69.0 

10.3 

Poisson's Ratio: 

v, 2 

0.3 

0.27 


0.3 

0.27 


0.3 

0.20 

Shear Moduli ( 10 9 Pa): 

G !2 

26.5 

7.20 

(j>2i 

26.5 

4.29 

G 3I 

26.5 

7.20 

Piezoelectric Charge Constant (10 12 

^31 

m/V) : 

-154. 

— 

Electric Permittivity ( 10‘ 9 f/m): 

£ 11 

15.05 

— 

£>22 

15.05 

— 


15.05 

— 

Thermal Expansion (10' 6 m/m°C): 

“n 

1.2 

-0.9 

a 22 

1.2 

27.0 

Reference Temperature, T c (°C): 

20.0 

20.0 
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Figure 12: Comparison of centerline sensory deflections of a clamped [0 8 /p] 
Carbon/Epoxy plate with a piezoceramic layer under a uniform 50 °C load. 
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Figure 13: Comparison of centerline sensory electric potentials of a clamped [0 8 /p] 
Carbon/Epoxy plate with a piezoceramic layer under a uniform 50°C load. 
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Figure 14: Comparison of centerline deflections of a clamped [0 8 /p] Carbon/Epoxy 
plate with a piezoceramic layer under a uniform 50 °C load with applied 
active voltages of 0, 250, and 500 V. 
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6.4 Additional Numerical Studies 


Results from numerical studies are presented to demonstrate additional capabilities of the 
current formulation to model the behavior of piezoelectric composite structures. The capabilities 
investigated include: incorporating pyroelectric effects and temperature dependent material 
properties into the analysis and evaluating the significance of these effects on the thermal response; 
demonstrating capabilities to achieve selective shape control of thermally induced bending and 
twisting deformations; assessing the improved accuracy of layerwise models in predicting thermal 
stresses; and determining the influence of curvature on the sensory and active responses. 

6.4.1 Pyroelectric Effects 

Pyroelectric effects arise in piezoelectric material due to temperature dependent changes in 
the electrical properties. This behavior is analogous to the change in mechanical strains observed in 
elastic materials due to thermal expansion. In order to demonstrate the effect of incorporating 
pyroelectric effects, the clamped plate problem shown in Figure 1 1 is re-examined incorporating 
pyroelectric effects. The problem details remain the same as before except for the consideration of 
non-zero pyroelectric constants in the piezoelectric material properties as listed in Table 6. The 
influence of the pyroelectric effects on the centerline deflection of the plate when the piezoelectric 
layer is operating in a sensory mode is shown in Figurel5 along with the corresponding result when 
the pyroelectric constant is neglected. The corresponding results for the sensory electric potential 
which develops is shown in Figure 16. Both figures demonstrate that pyroelectric effects can have 
a significant impact on both the displacements and electric potentials of piezoelectric composite 
structures in thermal environments. 


NASA/TM — 200 1-210892 


66 



Table 6: Material properties of piezoceramic and Carbon/Epoxy. 


Property 

Piezoceramic 

Carbon/Epoxy 
(Daniel and Ishai, 1994) 

Elastic Moduli ( 10 9 Pa): 

E„ 

69.0 

142.0 

E,, 

69.0 

10.3 

E„ 

69.0 

10.3 

Poisson's Ratio: 



Up 

0.3 

0.27 

"O73 

0.3 

0.27 

v 3 i 

0.3 

0.20 

Shear Moduli (10 9 Pa): 

G n 

26.5 

7.20 


26.5 

4.29 

G 3] 

26.5 

7.20 

Piezoelectric Charge Constant (10 12 m/V) : 


d 3 . 

-154. 

— 

Electric Permittivity (10 9 

f/m): 


£ n 

15.05 

— 

677 

15.05 

— 


15.05 

— 

Thermal Expansion (10 6 

m/m°C): 


an 

1.2 

-0.9 

^22 

1.2 

27.0 

Pyroelectric Constant (10 

3 C/m 2 °C): 


P 3 

-0.2 

— 

Reference Temperature, T 0 (°C): 20.0 

20.0 
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Figure 15: Pyroelectric effects on the centerline sensory deflections of a 
clamped [0 8 /p] Carbon/Epoxy plate with an attached piezoceramic layer 
under a uniform 50 °C load. 
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Figure 16: Pyroelectric effects on the centerline sensory electric potentials of a 
clamped [0 8 /p] Carbon/Epoxy plate with an attached piezoceramic layer 
under a uniform 50°C load. 
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6.4.2 Temperature Dependent Material Properties 


The various elastic, dielectric, piezoelectric, and pyroelectric material properties of 
piezoelectric materials are influenced differently by temperature variations. In addition, piezoelectric 
materials also possess a characteristic limiting temperature, called the Curie temperature, beyond 
which the material loses its piezoelectric properties. Thus, in typical applications the operating 
temperature of piezoelectric materials must be considerably lower than the Curie temperature. In the 
current study, the temperature dependent variations of an AS4/350 1 -6 Carbon/Epoxy composite and 
a PZT-5A piezoceramic were implemented using a piecewise linear representation in the constitutive 
equations (Equations 23-26). The temperature dependence of the AS4/3501-6 moduli and a 22 are 
depicted in Figure 17, while the variation of the PZT-5A properties d 3l and e 33 are shown in Figure 
18, while a n , and p-, are shown in Figure 19. All properties are nondimensionalized with their 
respective room temperature values shown in Table 7. The simply supported plate with discrete 
piezoelectric patches examined in Figure 9 is used to demonstrate the effects of temperature 
dependent material properties. Linear thermal gradients of different values are applied through-the- 
thickness of the plate with the bottom surface (z/h =-0.5) fixed at 20°C and the upper surface (z/h 
= 0.5) varying from 45° to 170°C. 

The first case utilizes a combined active/sensory configuration of the piezoceramic patches, 
in which the upper patches are utilized as sensors and the lower patches are used as actuators with 
zero volts applied. All the piezoceramic patches are grounded on the surface in contact with the 
carbon/epoxy plate. The resulting thermally induced deflection at the center of the plate is shown in 
Figure 20 for different applied thermal gradients. The deflections are nondimensionalized with 
respect to the total laminate thickness (H = 1 .0 1 mm). The results show that when material properties 
are assumed to remain constant with temperature, a linear variation in deflection occurs. In contrast, 
a nonlinear variation in the deflection is observed when the temperature dependent properties are 
incorporated. Larger discrepancies occur between the two cases as larger thermal gradients are 
applied. The corresponding sensory voltages at the central piezoceramic patch on the upper surface 
of the plate is shown in Figure 21. As with the deflections, the assumption of constant material 
properties leads to a linear variation of the voltages, while nonlinearities are introduced when 
temperature dependent properties are modeled. Closer agreement is observed in this case. 

A fully active configuration of the piezoceramic patches is used to investigate the effect of 
temperature dependent material properties on thermal shape control. Equal voltages are applied to 
all the piezoceramics. The applied voltages required to equally minimize the deflection at the center 
of the plate to zero for different applied voltages are shown in Figure 22. Once again, a nonlinear 
variation is predicted by the temperature dependent case and an increasing discrepancy is observed 
as the thermal gradient increases. This discrepancy corresponds to the trend observed for the center 
deflection of the plate (Figure 20). Since larger deflections occur for the temperature dependent case, 
this in part justifies the larger applied voltages required to minimize the deflections. 
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Figure 17: Temperature dependence of AS4/3501-6 Carbon/Epoxy 
properties (Daniel and Ishai, 1 994). 
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Figure 18: Temperature dependence of PZT-5A properties d 31 and e 33 . 
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Figure 19: Temperature dependence of PZT-5A properties a n and p 3 . 
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Table 7: Material properties of piezoceramic (PZT-5A) and Carbon/Epoxy 

(AS4/3501-6) composite. 


Piezoceramic 

Carbon/Epoxy 
(Daniel and Ishai, 1994) 

Elastic Moduli (10 9 Pa): 

E„ 

69.0 

142.0 

E-,-, 

69.0 

10.3 

e 33 

53.0 

10.3 

Poisson's Ratio: 

«12 

0.31 

0.27 

*>23 

0.44 

0.20 


0.38 

0.02 

Shear Moduli (10 9 Pa): 

0,2 

23.3 

7.2 

G 23 

21.1 

4.29 

0 3 l 

21.1 

7.2 

Thermal Expansion (10 6 m/m°C): 

a,. 

1.2 

-0.9 

a 22 

1.2 

27.0 

Piezoelectric Charge Constant (10 12 m/V) : 

d 31 

-154. 

— 

Electric Permittivity ( 10 9 f/m): 

£ 33 

15.05 

— 

Pyroelectric Constant (10 3 C/m 2 °C): 

P3 

-2.0 

— 

Reference Temperature, T 0 , (°C): 

20. 

20. 

Curie Temperature, T c , (°C): 

365. 

— 
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Figure 20: Influence of temperature dependent material properties on the deflections 
at the center of a [0/±45] s Carbon/Epoxy plate with 30 attached piezoceramic 
patches in an active/sensory mode with varying applied thermal gradients 
through-the-thickness (25 °C to 150°C) . 
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Figure 21: Influence of temperature dependent material properties on the electric 
potentials at the center of a [0/±45] s Carbon/Epoxy plate with 30 attached 
piezoceramic patches in an active/sensory mode with varying applied 
thermal gradients through-the-thickness (25 °C to 150°C) . 
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Figure 22: Influence of temperature dependent material properties on the active 
voltages to minimize deflections at the center of a [0/±45] s Carbon/Epoxy plate 
with 30 attached piezoceramic patches in an active/active mode with varying 
applied thermal gradients through-the-thickness (25 °C to 150 °C) . 
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6.4.3 Selective Shape Control of Bending and Twisting 


The potential of piezoelectric actuators to selectively minimize thermally induced 
deformations is demonstrated on a clamped [45 3 /-45,] plate with attached piezoelectric patches. The 
geometry and material properties remain the same as in Figure 9 and Table 4, respectively. The plate 
is clamped along x=0 and a 30°C thermal gradient is applied through-the-thickness with the 
piezoelectric patches operating as actuators. The temperature dependence of material properties is 
neglected in this study. The resulting initial deformation is shown in Figure 23 and shows a 
combined bending and twisting deformation induced by the thermal gradient. The bending 
deformation can be selectively minimized by applying the same voltages to the top and bottom 
piezoceramic patches. The resulting deformation when 40 V is applied is shown in Figure 24. The 
twisting deformation remains unaffected and is now more apparent due to the removal of the 
bending. In a similar manner, the twisting deformation can be minimized without affecting the 
bending behavior. This is accomplished by applying voltages of opposite polarities on the top and 
bottom surface piezoceramic patches. The deformation produced by the application of 145 V on the 
top piezoceramic patches and - 145 V on the bottom piezoceramic patches is shown in Figure 25. The 
bending deformation is more apparent due to the removal of the twisting. In order to minimize both 
the bending and twisting deformations simultaneously, the previous two cases can be superimposed. 
Thus, by applying 185 V on the top actuators and -105 V on the bottom actuators, all deformations 
in the plate can be eliminated as shown in Figure 26. 
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Figure 23: Initial bending and twisting deformation of a clamped [45 3 /45. 3 ] 
AS4/3501 plate.with 30 attached discrete G-1195 piezoelectric patches 
under a 30 °C thermal gradient through-the-thickness. 
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Figure 24: Selective shape control of the bending deformation of a clamped [45 3 /45. 3 ] 
AS4/3501 plate. with 30 attached discrete G-l 195 piezoelectric patches 
under a 30 °C thermal gradient through-the-thickness by applying 
40 V to all piezoceramic patches. 
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Figure 25: Selective shape control of the twisting deformation of a clamped [45/45.-,] 
AS4/3501 plate. with 30 attached discrete G-l 195 piezoelectric patches 
under a 30 °C thermal gradient through-the-thickness by applying 
145 V to top patches and -145 V to bottom patches. 


NASA/TM — 200 1 -2 1 0892 


81 



15 -r- 

10 + 



Figure 26: Shape control of both bending and twisting deformation of a clamped 
[45 3 /45. 3 ] AS4/3501 plate.with 30 attached discrete G-l 195 piezoelectric patches 
under a 30 °C thermal gradient through-the-thickness by applying 
185 V to top patches and -105 V to bottom patches. 
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6.4.4 Thermal Stresses 


The increased accuracy gained from implementing a layerwise laminate theory is 
demonstrated by examining the thermal stresses in a cantilevered [0 8 /p] beam. The beam is depicted 
in Figure 27 and has a length (L) of 25.4 cm and a width (b) of 2.54 cm. Two thicknesses (H) are 
examined; a thin beam (H = 0. 1 143 cm) with an aspect ratio (L/H) of -200 and a thick beam (H = 
11.43 cm) with an aspect ratio (L/H) of -2. The beam is composed of eight unidirectional E- 
Glass/Epoxy plies and one piezoelectric layer. The material properties are listed in Table 8. The 
beam element is used to conduct the analysis using 30 elements along the length of the beam, with 
one discrete layer modeled through-the-thickness for each piezoelectric and E-Glass/Epoxy ply. A 
uniform temperature rise of 100°C is applied with the piezoelectric layer grounded (i.e. 0 V applied). 
The predicted normal stress, o,„ near the middle of the beam (x/L = 0.5) is show'n in Figure 28 for 
both the thin and thick beams, along with the corresponding results predicted by classical laminate 
theory (CLT). The normal stress is nondimensionalized using the equivalent laminate modulus, E^, , 
= 42.2 10 9 Pa. The CLT predicts the same stress regardless of the thickness of the beam. At this 
location of the beam, there is good agreement between the current formulation and CLT regardless 
of the beam thickness. The predicted normal stress, c M , at the free end of the beam (x/L = 1.0) is 
shown in Figure 29 and shows that CLT is only accurate for the thin beam. The predicted shear 
stress, a 13 , at the free end of the beam (x/L = 1.0) is shown in Figure 30. The shear stress is 
nondimensionalized using the equivalent laminate modulus, G L13 = 6.3 10 9 Pa. The results show that 
the CLT assumption of neglecting shear stresses is only accurate for the thin beam case. These 
results help demonstrate the improved accuracy in stress predictions achieved by using the layerwise 
laminate theory, especially for thick laminates. 
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Figure 27: Cantilevered [0 p /p] E-Glass/Epoxy and 
piezoelectric beam 
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Table 8: Material properties of piezoelectric and E-Glass/Epoxy. 



Piezoelectric 

E-Glass/Epoxy 
(Daniel and Ishai, 1994 ) 


Elastic Moduli (10 9 Pa): 

E„ 

68.0 

39.0 


E33 

68.0 

8.6 


Poisson’s Ratio: 


0.30 

0.28 


Shear Modulus (10 9 Pa): 

^13 

26.2 

3.8 


Thermal Expansion (10 6 m/m°C): 

On 

3.8 

7.0 


<*33 

3.8 

21.0 


Piezoelectric Charge Constant (10 

d 31 

12 m/V) : 

-125. 

— 


Electric Permittivity (10 9 f/m): 

£ 33 

11.06 

— 


Pyroelectric Constant ( 1 O' 3 C/m 2 0 
P3 

C): 

-2.5 

— 


Reference Temperature, T 0 , (°C): 

20. 

20. 
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Figure 28: Comparison of normal stresses at the middle (x/L=0.5) of a [0 8 /p] 
E-Glass/Epoxy and piezoelectric beam under a uniform 100°C load for two 
aspect ratios (thin=200 and thick=2). 
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Figure 29: Comparison of normal stresses at the free end (x/L=1.0) of a [0 8 /p] 
E-Glass/Epoxy and piezoelectric beam under a uniform 100°C load for two 
aspect ratios (thin=200 and thick=2). 
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Figure 30: Comparison of shear stresses at the free end (x/L=1.0) of a [0 8 /p] 
E-Glass/Epoxy and piezoelectric beam under a uniform 100°C load for two 
aspect ratios (thin=200 and thick=2). 
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6.4.5 Influence of Curvature on Active and Sensory Response 

The problem examined consists of a [0,/p] (where 0° is along the \ axis) circular cylinder 
with a radius (R) of 0.76 m, a length (L) of 1 .524 m, and a thickness (h) of 7.62 mm as depicted in 
Figure 31. The cylinder is composed of nine unidirectional Carbon/Epoxy plies and one 
piezoceramic layer. The material properties are listed in Table 6. Due to the symmetry of the 
problem, only Vs of the cylinder is modeled using a 10x10 mesh with two discrete layers (one for the 
Carbon/Epoxy and one for the piezoceramic). A clamped-clamped configuration in which both ends 
of the cylinder (z/L = 0 and z/L = 1 ) are fixed is examined. The cylinder is subjected to three 
different types of thermal loads (a uniform thermal load and two types of sinusoidally varying 
temperatures along the hoop direction) as shown in Figure 32. Results from the current shell element 
are presented to investigate the influence of the curved structure on the sensory and active response 
under the different applied thermal loads. All displacements (w) and electric potentials (<p) presented 
in this section are along the hoop direction at z/L = 0.5 and r/R =1.0 of the cylinder. 

The influence of pyroelectric effects on the displacements and sensory electric potentials of 
the cylinder under the three different applied thermal loads are presented in Figures 33-38. Each 
figure contains results for two cases to show the effects of incorporating the pyroelectric constant 
in the analysis (the p 3 =0 case neglects pyroelectric effects). Figures 33 and 34 show, respectively, 
the displacement and sensory electric potential induced from applying a uniform thermal load ol 
50°C. A uniform deflection of the cylinder is achieved, which maintains the original circular shape. 
The sensory electric potential also shows a uniform response that corresponds to the displacement. 
The incorporation of pyroelectric effects produces larger displacements and electric potentials. The 
results from applying a cosine varying temperature are shown in Figures 35 and 36. The deflection 
of the cylinder now shows a sinusoidal pattern that translates into either an oval-shaped or figure 
eight shaped circumferential displacement, depending on whether pyroelectric effects are modeled. 
The electric potential only displays a sinusoidal pattern when pyroelectric effects are included and 
becomes almost zero when the pyroelectric effect is neglected. Figures 37 and 38 show the results 
of applying a double cosine varying thermal load. A sinusoidally varying displacement and electric 
potential pattern is produced which corresponds to an oval-shaped deflection of the cylinder. Once 
again, incorporating the pyroelectric effect leads to increased displacements and electric potentials. 

The active response of the cylinder under the three applied thermal loads is examined in 
Figures 39-41. Each figure shows the displacement of the cylinder for three configurations of the 
piezoelectric layers: a sensory mode, a grounded configuration (zero electric potential applied), and 
an active mode. The magnitude of the applied electric potential remains the same for all three cases, 
although the form of the electric potential varies to correspond to the type of thermal load. Figure 
39 shows the deflection of the cylinder under a uniform thermal load (50°C). Application of 
increasing uniform applied electric potentials produces a noticeable decrease in the thermally 
induced sensory deflection. Figure 40 depicts the displacement of the cylinder under a cosine 
varying thermal load and shows that although some changes in the overall shape of the cylinder can 
be obtained by applying active electric potentials, no significant overall reduction of the thermal 
deflection is achieved. Figure 41 illustrates the deflections produced under a double cosine variation 
of the thermal load. For this case, the double cosine applied electric potential effectively minimizes 
the thermal deflection. 
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Figure 31 : Geometry of a [0</p] Carbon/Epoxy circular cylinder with an 
attached piezoceramic layer. 
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Figure 32: Different types of applied thermal loads applied to the cylinder: 
a) 50°C uniform thermal load; b) 50°C cosine varying load; and 
c) 25 °C double cosine varying load. 
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Figure 33: Influence of pyroelectric effects on the displacements of a [0</p] 
Carbon/Epoxy circular cylinder with a piezoceramic layer under 
a 50 °C uniform thermal load. 
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Figure 34: Influence of pyroelectric effects on the sensory electric potentials of a [0,/p] 
Carbon/Epoxy circular cylinder with a piezoceramic layer under 
a 50 °C uniform thermal load. 
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Figure 35: Influence of pyroelectric effects on the displacements of a [0</p] 
Carbon/Epoxy circular cylinder with a piezoceramic layer under 
a 50°C cosine varying thermal load. 
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Figure 36: Influence of pyroelectric effects on the sensory electric potentials of a [0,/p] 
Carbon/Epoxy circular cylinder with a piezoceramic layer under 
a 50 °C cosine varying thermal load. 
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Figure 37: Influence of pyroelectric effects on the displacements of a [0</p] 
Carbon/Epoxy circular cylinder with a piezoceramic layer under 
a 25 °C double cosine varying thermal load. 
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Figure 38: Influence of pyroelectric effects on the sensory electric potentials of a [(yp] 
Carbon/Epoxy circular cylinder with a piezoceramic layer under 
a 25 °C double cosine varying thermal load. 
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Figure 39: Comparison of displacements of a [0^/p] Carbon/Epoxy circular cylinder 
with a piezoceramic layer in sensory, grounded, and active modes 
under a 50 °C uniform thermal load. 
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Figure 40: Comparison of displacements of a [0</p] Carbon/Epoxy circular cylinder 
with a piezoceramic layer in sensory, grounded, and active modes 
under a 50 °C cosine varying thermal load. 
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Figure 41: Comparison of displacements of a [0g/p] Carbon/Epoxy circular cylinder 
with a piezoceramic layer in sensory, grounded, and active modes 
under a 25 °C double cosine varying thermal load. 
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CHAPTER 7 

CONCLUSION AND SUGGESTIONS FOR FUTURE WORK 


New analytical formulations and mechanics are developed to characterize the coupled 
mechanical, electrical, and thermal behavior of piezoelectric composite materials. The coupled 
thermopiezoelectric behavior is captured at the material level through the constitutive equations and 
the displacements, electric potential, and temperature are introduced as state variables in the 
formulation. The equations of motion are developed in both cartesian and curvilinear coordinates 
and a layerwise laminate theory is incorporated to provide more accurate predictions of interlaminar 
and intralaminar effects, especially for thick laminates. Thermal effects arising from thermal 
expansion mismatch, pyroelectric effects, and temperature dependent material properties are 
incorporated in the formulation. Corresponding finite element equations are developed and 
implemented into stand alone computer programs for beam, plate, and shell elements. The accuracy 
of the formulation is verified with comparisons from published experimental data and other 
analytical methods. Additional numerical studies are also conducted to demonstrate additional 
capabilities of the formulation to predict the sensory and active behavior of piezoelectric composite 
materials. The verification and numerical studies have helped to demonstrate the effectiveness and 
accuracy of the developed mechanics and finite elements. However, future experimental studies will 
help to complement this work by contributing to additional verification studies. 

7.1 Suggested Plan of Future Experimental Studies 

In order to verify the accuracy of the developed analytical formulation in more detail, the 
following plan of experimental studies are outlined. As illustrated in the previous sections, previous 
researchers have conducted numerous room temperature static experimental studies along with some 
limited high temperature studies. However, dynamic studies, and in particular high temperature 
dynamic experiments, have received less attention. Thus, the development of an experimental setup 
for conducting high temperature dynamic studies will help to more fully characterize the behavior 
of piezoelectric composite materials and determine the feasibility of applying these materials in 
extreme temperature environments. 

A generic dynamic experimental setup is shown in Figure 39. A beam specimen with an 
attached piezoelectric patch will be excited by an electromagnetic shaker and the corresponding 
response of the beam measured by an accelerometer. For high temperature tests, the beam will be 
placed inside an oven or furnace. The time varying responses of the beam and piezoelectric patch 
are passed through various test measurement devices including signal conditioners, amplifiers, and 
oscilloscopes to increase the signal strength and gain, as well as to monitor and measure the dynamic 
waveforms. The time signals are then passed on to a PC based dynamic data acquisition system 
which contains hardware to acquire the signals and software to calculate the frequency response of 
the data. The resulting output will include the natural frequencies of the beam, the mode shapes, and 
the sensor output. 

A series of tests will be performed to investigate the individual sensory and active dynamic 
responses of piezoelectric composite materials at high temperature. The first test configuration will 
consist of a beam with a single piezoelectric patch as shown in Figure 40. A series of tests will be 
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performed as listed in Table 9 to determine the sensory and active responses individually. A room 
temperature baseline test will be performed to determine the natural frequencies of the beam with 
the piezoelectric patch grounded. The piezoelectric patch will then be tested in a sensory mode to 
demonstrate the capability to monitor modal frequencies through measurement of the sensory 
voltages, and to determine the effect of the sensory mode on the natural frequencies. Next, the 
piezoelectric patch will be configured as an actuator and electric voltages applied to measure the 
capability to actively control the vibration of the beam. These experiments will then be repeated at 
higher temperature to assess thermal effects on the sensory and active responses. A series of different 
beam materials from metals to composites can be studied, as well as different types of piezoelectric 
materials from piezopolymers to piezoceramics, in order to determine the most effective sensor and 
actuator materials. In addition, the size and location of the piezoelectric patch can be varied to assess 
the influence of coverage area and location on sensing and actuation. 

The next test configuration will consist of a beam with two piezoelectric patches, as shown 
in Figure 41, to investigate combined sensory/active and active/active modes of operation. A series 
of tests will be performed as shown in Table 10 to determine capabilities for combined sensory and 
actuator operations. The baseline case will consist of a room temperature test with both piezoelectric 
patches grounded to determine the natural frequencies of the beam. A combined sensory/active 
configuration will then be studied, in which various voltages will be applied to the actuator patch 
to reduce the natural frequencies, while the sensor patch monitors the response. The next case will 
involve both patches operating as actuators to determine whether multiple actuators can more 
effectively reduce the beam vibrations than a single actuator. These experiments will then be 
repeated at higher temperature to determine thermal effects on the combined sensory/active and 
active/active configurations. Additional cases with different beam and piezoelectric materials, as 
well as different sized patches and locations will also be studied to determine the optimal 
combinations for combined sensory/active and active/active operations. 

Additional studies will also be conducted on a flat plate with multiple piezoelectric patches 
(Figure 42) and a curved plate with multiple piezoelectric patches (Figure 43). The flat plate will 
provide additional information on two-dimensional geometry effects on sensing and actuation, as 
well as to determine capabilities to utilize patches at specific locations to sense and actuate certain 
natural frequencies. A similar study will also be performed on a cylindrical shell to determine the 
effect of curvature on actuation and sensing. As in previous problems, the baseline case will involve 
determining the room temperature natural frequencies with the piezoelectric patches inoperative. A 
variety of different sensor and actuator configurations will be examined at room and high 
temperatures, for different beam and piezoelectric materials. 

Throughout the experimental studies, the developed finite element codes will be used to 
provide design guidelines to reduce the number of experiments. Even for the simple beam 
configurations studied in Figures 40 and 4 1 , there a large number of different experiments required 
to determine the influence of different materials, patch sizes, and locations. As more complicated 
geometries are considered in Figures 42 and 43, the number of experimental tests will increase 
substantially. Thus, in order to keep experimental testing time and costs reasonable, a combined 
analytical and experimental approach is essential. 
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Figure 39: Dynamic experimental setup. 
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Testing Configuration 



Figure 40: High temperature test configuration for a cantilevered beam 
with one piezoelectric patch to examine sensory and 
active modes individually. 
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Table 9: Series of tests for a cantilevered beam with one piezoelectric patch. 




Beam and Piezoelectric Materials 



Aluminum Beam 

[0 8 ] Graphite/Epoxy 
Beam 

tom], 

Graphite/Epoxy Beam 

Temp. 

Piezo 

Mode 

PZT 

Patch 

PVDF 

Patch 

PZT 

Patch 

PVDF 

Patch 

PZT 

Patch 

PVDF 

Patch 

Room 

Grounded 







Room 

Sensory 







Room 

Active 







Elevated 

Grounded 







Elevated 

Sensory 







Elevated 

Active 
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Testing Configuration 



Figure 4 1 : High temperature test configuration for a cantilevered beam 
with two piezoelectric patches to examine combined sensory/active 
and active/active modes of operation. 


N AS A/TM— 200 1 -2 1 0892 


106 





Table 10: Series of tests for a cantilevered beam with two piezoelectric patches. 




Beam and Piezoelectric Materials 



Aluminum Beam 

[0 g ] Graphite/Epoxy 
Beam 

[0 2 /90 2 

| s Graphite 
Beam 

:/Epoxy 

Temp. 

Piezo. 

Mode 

PZT/ 

PZT 

PVDF/ 

PVDF 

PZT/ 

PVDF 

PZT/ 

PZT 

PVDF/ 

PVDF 

PZT/ 

PVDF 

PZT/ 

PZT 

PVDF/ 

PVDF 

PZT/ 

PVDF 

Room 

Grounded 










Room 

Active/ 

Sensory 










Room 

Active/ 

Active 










High 

Grounded 










High 

Active/ 

Sensory 










High 

Active/ 

Active 
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Piezoelectric patch 



Figure 42: Test configuration of a flat plate with attached piezoelectric 
patches to investigate two-dimensional effects on combined 
sensory and active operations. 
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Figure 43: Test configuration of a curved cylinder with attached piezoelectric 
patches to investigate curvature effects on combined 
sensory and active operations 
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